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O '. Abstract 

Q 

^+ ■ 

The emission rate of photons from a hot, weakly coupled ultrarelativistic 
plasma is analyzed. Leading-log results, reflecting the sensitivity of the emis- 
^ ■ sion rate to scattering events with momentum transfers from gT to T, have 

'nI" . previously been obtained. But a complete leading-order treatment requires 

, including collinearly enhanced, inelastic processes such as bremsstrahlung. 

0^ \ These inelastic processes receive 0(1) modifications from multiple scattering 

during the photon emission process, which limits the coherence length of the 
emitted radiation (the Landau-Pomeranchuk-Migdal effect). We perform a 
(■H , diagrammatic analysis to identify, and sum, all leading-order contributions. 

J^, We find that the leading-order photon emission rate is not sensitive to non- 

Qh! perturbative g^T scale dynamics. We derive an integral equation for the pho- 

ton emission rate which is very similar to the result of Migdal in his original 
discussion of the LPM effect. The accurate solution of this integral equation 
for specific theories of interest will be reported in a companion paper. 
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I. INTRODUCTION AND RESULTS 

The recent commissioning of RHIC, and anticipated heavy ion experiments at the LHC, 
have stimulated work on photon emission in both equilibrium and nonequilibrium relativis- 
tic plasmas [p!|-[TT[|. At sufficiently high temperatures or energy densities, QCD interactions 



become weak, and perturbation theory can be applicable for observables which are predom- 
inantly sensitive to the dynamics of typical excitations in the plasma, those with momenta 
comparable to the temperature. Whether perturbation theory is really useful for the energy 
densities achievable at RHIC and the LHC is an open question, but thermal field theory also 
has important applications in early universe cosmology (such as baryogenesis, leptogenesis, 
and magnetogenesis) . From a purely theoretical perspective, hot, weakly coupled, relativis- 
tic plasmas provide an instructive and interesting domain with rich dynamics dependent on 
the interplay of a wide variety of quantum and statistical effects, in which one can, with 
effort, perform controlled calculations for physically interesting quantities. 



A. Background 



The basic perturbative tools for studying a relativistic plasma near equilibrium were 
developed long ago [Q. In the last decade or so, there has been substantial progress. It 
is well understood today that there are three distinct length scales on which a hot weakly 
coupled relativistic plasma exhibits different characteristic behavior. 

The first scale is the "hard" scale, corresponding to wavenumbers (or momenta) A; ~ T. 
This is the characteristic scale of momentum or energy for the vast majority of the excitations 
(quarks, gluons, ...) comprising the plasma. The contributions of such excitations dominate 
bulk thermodynamic properties such as the total energy and momentum density of the 
plasma, as well as conserved charge susceptibilities (electric charge and various conserved or 
approximately conserved flavor charges, such as baryon number and isospin). The behavior 
of these excitations also controls transport properties such as viscosity, baryon number 
diffusion, and electrical conductivity |]T^|-|T^. Provided the gauge coupling^ is weak, (7 -C 1, 
excitations with wave number ~ T propagate as almost free quasi-particles moving at 
essentially the speed of light, engaging in occasional small angle scatterings [deflections by 
an angle of 0{g) or less] with a mean free time of order ^ ^/T, and even more 

occasional large angle scatterings [with 0(1) deflection] on times scales of order 1/g^T. 
These parametric estimates ignore factors of lng~^ which are present in both mean free 



times 14-17 



The second scale is the "soft" (or "semi-hard") scale, k ~ gT. For excitations of this 
wave number, plasma effects such as Debye screening, plasma oscillations, and Landau 
damping, become of 0(1) importance. These arise because of coherent interactions of the 
gT scale fields with the thermal bath of harder, 0(T) scale excitations. Dispersion relations, 
and non-Abelian gauge interactions, receive 0(1) corrections which must be appropriately 
resummed into the propagators and vertices [p!8| -|20[|. This resummation is commonly referred 
to as "hard-thermal-loop" (HTL) perturbation theory. After this selective resummation, 
the effective interactions amongst the gT scale degrees of freedom are weak, and hence 
perturbation theory remains applicable for gT scale dynamics. However, perturbation theory 
at the gT scale is only an expansion in g, not g^. The primary physical importance of gT 
excitations comes from the role they play as exchanged particles in scattering processes. The 
"transport" (or large-angle) scattering rate is dominated by exchange momenta in the range 
from T down to gT, while the total scattering rate is dominated by exchange momenta in 
the range from gT down to g^T. 

The third scale is the "ultrasoft" or non-perturbative scale, k ~ g'^T. For gauge bosons 
(or scalars, if the temperature is sufficiently close to a phase transition temperature), the 
occupation number of individual modes at this scale is so large, 0(1/(7^), that the mutual 
interactions between g'^T scale degrees of freedom are truly strongly coupled; perturbative 
treatments break down. This scale is important for problems which require non-perturbative 



"'^Here, and throughout our discussion, the gauge couphng is to be understood as the relevant couphng 
defined at the scale of the temperature, g = g{T). For hot QCD, g means the strong gauge coupling gs, 
while for hot QED, g should be understood as the electric charge e. Our treatment is applicable to both 
cases, although a few of the following comments concerning g^T scale dynamics assume that the gauge field 
is non-Abelian. 
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FIG. 1. The two processes of interest, bremsstrahlung and inelastic pair annihilation. In each case, the 
emerging photon is hard, with energy ~ T, but is nearly coUinear with the quark from which it is radiated. 
One or more interactions which exchange momentum with other excitations in the plasma are required for 
these processes to occur; the exchange momenta are all soft, with energy and momentum ~ gT. Time should 
be viewed as running from left to right. 



physics, such as baryon number violation in the hot standard model above the electroweak 



phase transition 21 



One would expect that questions involving the long time behavior of sensible observ- 
ables, dominated by A; ~ T excitations, should only depend on T and gT scale physics, at 
least at leading order in powers of g, and should therefore be perturbatively computable.0 
Numerous physically interesting questions fall into this category. Some examples are the 
hard photon emission rate from a thermalized quark-gluon plasma, the energy loss of a very 
hard parton (thermal "jet quenching"), and transport coefficients such as shear and bulk 
viscosity, electrical conductivity, and various diffusion constants (baryon number, isospin, 
etc.). In general, the question of how fast a nonequilibrium plasma relaxes to equilibrium, 
and which degrees of freedom have the slowest approach to equilibrium, depends on the 
dynamics of typical hard excitations. 

Surprisingly, not one of the physical observables just mentioned has been correctly com- 
puted, in a hot gauge theory, beyond leading order in the logarithm of the coupling, \og{l/g) 
(that is, neglecting corrections suppressed only by l/\ng~^). The chief obstruction to such 
an evaluation is the fact that each of these problems depends, at leading order in ^f, on the 
processes depicted in Fig. |I], namely bremsstrahlung and inelastic pair annihilation, with 
soft exchange and nearly collinear external states. These processes contribute to trans- 
port mean free paths and to the photo-emission rate at leading order, despite being 2^3 
particle processes, because a combination of soft and collinear enhancements make them 
occur at an 0{g'^T'^) rate per spacetime volume [or, for photo-emission from a quark gluon 
plasma, O^e^g^T'^)], which is the same as 2 ^ 2 scattering processes. The importance of 
these processes have been emphasized recently by Aurenche et al. 0, who demonstrated 
that they contribute at leading order to photon emission from the quark-gluon plasma. This 
was an important realization, but the quantitative analysis was incomplete (as recognized 
by the same authors in 0,|]), because it ignored a suppression of the emission rate due 
to the Landau- Pomeranchuk-Migdal (LPM) effect. That is, one internal propagator in the 



^The total scattering rate for a hard excitation, with its logarithmic sensitivity to g^T scale dynamics, is 
not really a meaningful observable, for roughly the same reason that the strictly elastic scattering amplitude 
for electrons in QED is not meaningful. 
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diagram is very nearly on-shell and receives 0(1) corrections from the quasiparticle width. 
This reflects sensitivity to additional scatterings which occur during the photon emission 
process. However, merely including the appropriate width on the intermediate propagator 
is not sufficient; a consistent treatment requires a more elaborate and detailed analysis. The 
purpose of this paper is to provide this analysis, and to derive an integral equation, similar in 
form to a linearized kinetic equation, whose solution determines the rate of bremsstrahlung 
and inelastic pair annihilation processes, to leading order in g. The accurate solution of this 



integral equation for specific theories of interest will be reported in a companion paper [22 



At times in our discussion we will refer to the relevant charge carriers as "quarks" and 
the exchanged gauge bosons as "gluons" , but it should be emphasized that our results apply 
equally well to hot QED. We assume that the relevant interactions in the plasma come from 
a weakly coupled simple gauge group, either Abelian or non-Abelian. We also assume that 
the emitted photon has negligible interactions with the plasma subsequent to its emission.^ 
Although it is also physically quite interesting, we will not consider the case where the hard 
coUinear emitted particle is a gluon instead of a photon, but leave this for future work. 



B. Soft exchange with collinear emission 

Let us first indicate why the processes shown in Fig. |1| are important, by sketching the 
evaluation of the bremsstrahlung rate; estimating the pair annihilation rate is completely 
analogous. Our argument roughly follows that of Aurenche et al. . Assume, for simplicity, 
that the photon momentum is hard, k T. The analysis is easiest for scalar quarks, although 
the conclusions we will find also hold for real fermionic quarks. Labeling the momenta as 
shown in the figure,^ and requiring that the outgoing states be on-shell fixes K'^ = and 
Pi = P2 = L\ = L\ = — ~ —g'^T'^, with the asymptotic thermal mass (squared) 
of the quarks. (Scalars have a thermal mass which is momentum-independent at one-loop 
order, while fermions acquire a momentum-dependent thermal dispersion relation. But for 
P ~ T, the momentum dependence is sub- leading in g and it is sufficient, at leading order, 
to just retain the momentum independent asymptotic thermal mass moo. This asymptotic 
mass is the same as the thermal mass for scalars provided the scalar quartic coupling is 
small compared to g'^.) 

If the exchange momentum Q is of order T, then the rate for this process is 0{e^g^), 
which is negligible compared to the 2^2 Compton scattering and annihilation processes 



•^Of course, the propagation of a photon really is affected by the medium it traverses; in particular, thermal 
corrections to the photon dispersion relation increase the energy of a photon of momentum ~ T by an 
amount of order e^T. In a hot QCD plasma, the photon coupling e is small compared to the strong coupling 
constant gs, so thermal corrections to the photon dispersion relation are negligible compared to the (non- 
electromagnetic) thermal effects we will be focusing on. For QED, instead of speaking of the photon emission 
rate, which ignores rescattering, one should consider the photon emissivity. Our analysis is also applicable in 
this case, although one would need to solve the final integral equation using a timelike photon 4-momentum 
suitably determined by its thermal dispersion relation. 

*As is customary in thermal field theory, we work in the plasma rest frame, write 4-momenta in capitals 
(Q), their spatial components in bold face (q), the magnitude of the three-momentum in lower case {q = |q|), 
and the time component in lower case with zero index {q^). 
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shown in Fig. whose rates are 0{e^g'^) up to logs p^j2^ . But for Q ~ gT there is an 
enhancement from the soft gluon propagator. The leading contribution will arise from the 
kinematic regime where g° ~ |q| ~ gT and Pi ■ Q ~ g'^T'^- The latter condition restricts 
the coUinearity of the outgoing photon and quark. Generically, if all components of Q are 
0{gT) then Pi - Q ^ gT"^; however this region turns out to be subdominant. (See footnote |^ 
below.) Because Q is small and the other 4-momenta are almost lightlike, the condition 
Q ■ Pi ^ g'^T'^i together with the requirement that p\ be at most 0{T), also implies that 
Q ■ P2, K ■ P2, and K ■ Pi are all 0{g^T'^), and restricts the spatial components of pi, p2, 
and q orthogonal to k to be 0{gT). 

The emission amplitude for a transverse photon with polarization vector e'^ is 



{P2 + K)^ + m'^ ' " ' {Pi - K)^ + m\ 

(1.1) 



where the two terms in the bracket come from the diagram where the photon is radiated 
from the outgoing P2 line (shown in Fig. |I]) and the corresponding diagram where the photon 
is radiated from the incoming Pi line. 

One may choose the photon polarization vectors to be purely spatial in the plasma frame 
and orthogonal to k. Therefore the quantities e ■ {Pi+P2±Q) appearing above are 0{gT). 
There is no special relation between the two terms which might cause a cancellation,^ so 
it is sufficient to consider the square of one term to estimate the rate for the process. 
There is an explicit e^g^ plus a g'^ phase space suppression (from the conditions that all 
components of q be small and that Q ■ Pi ~ g'^T'^)- The square of the soft gauge boson 
propagator G^etlQ) is of order l/Q"^ ~ 1/g'^T'^, and the square of the scalar propagator is 
1/[{P2+Kf + m^]2 ~ l/g^T\ Both (La+Li)^ and (Pi+Ps-iT)^ are 0(T). Finally there 
are the squares of the numerators, [(P1+P2+Q) -e]^ ~ Q^T^- Adding powers of g, we find the 
rate for the process is 0{e^g'^T^), which is the same order as the rate of the 2^2 processes 



^In the generic soft exchange regime where Q ■ Pi ^ gT^, the transverse (to k) components of Pi and 
P2 are order and equal [up to 0{gT)]. Hence, (Pi+P2±Q) • e ^ (P1+P2) • e g^^'^T. In this case, 

the two terms in Eq. (1.1) are equal and opposite at leading order and there is a cancellation, related to a 
cancellation addressed by Bodeker in |2^. This is why the more highly coUinear region with Q ■ Pi ^ g^T"^ 
is actually dominant. 
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shown in Fig. ^ Hence, the inelastic 2^3 processes of Fig. |] are of 0(1) importance in 
hard photon emission.^ 

Several points about the collinear emission, soft exchange regime which dominates these 
2^3 processes should be emphasized: 

1. The large contribution from a nearly on-shell quark propagator, 1 / {P2 + K)'^ ~ 
is essential. 

2. The large contribution from a gauge boson propagator with soft exchange momentum, 
1/Q2 ^ 1/^2^-2^ ig essential. 

3. The condition Pi • Q ~ g'^T'^, combined with near-lightlike on-shell incoming and 
outgoing momenta, forces P2 ■ K g'^T'^, which implies that the outgoing photon and 
quark are highly collinear. The opening angle between the outgoing quark and photon 
satisfies 1 — cos(6') ~ g"^, or 9 = 0{g). 

4. Although the above estimates used scalar charge carriers for simplicity, all of the 
conclusions are equally applicable to fermions. 

However, carrying out the calculation of the 2 3 process as we have just outlined 
is not sufficient to obtain the correct leading order emission rate. To see this, consider 
the first point above. The real part of the inverse scalar propagator appearing in our 
(sketch of a) calculation is 0{g'^T'^). But hard, charged particle self-energies also have 
0{g'^T^) imaginary parts. Therefore it appears necessary to include this imaginary part 
in the propagator. But the imaginary part characterizes the probability that the particle 
undergoes additional scattering between the time of the initial scattering shown in the 
diagram and the time of the photon emission. If it is necessary to include the imaginary 
part in an intermediate propagator, then it is presumably also important to include such 
scatterings explicitly, meaning that we must consider 3^4 processes, 4 5 processes, or 
more. 

Alternately, consider point 3 above. The collinearity condition requires that the com- 
ponents of the spatial vectors p2 and k which are transverse to pi be at most 0{gT). To 
resolve momenta with that precision, the wave packets for the outgoing states must have 
spatial extent of order 1/gT in the transverse direction. Since the opening angle is order g, 
the outgoing quark and photon must overlap in space over a time interval of order 1/g'^T. 
This is the formation region; the photon is generated from coherent interaction between the 
incoming and scattered waves over a region of this length, as roughly depicted in Fig. |^. 
However, the mean free path for soft scatterings is 0{l/g'^T), which is the same size. There- 
fore it will be necessary to consider corrections to the emission process due to scatterings 
occurring before the photon formation is complete. Such scatterings can potentially disrupt 
the photon emission process; this is called the Landau-Pomeranchuk-Migdal effect [p^30|. 



Because the emission time is 0{l/g^T), only scattering processes with a mean free path of 
order 1/g'^T will be important. Therefore we can neglect scattering processes with exchange 



^ Since the incoming and outgoing quark momenta, Pi and P27 differ by 0{T), tlie analogous pure QCD 
processes involving hard gluon emission are also of 0{\) importance in transport coefficients, as recently 
emphasized in P6| . 
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FIG. 3. A cartoon of the real-space appearance of scattering and photon emission. The photon emission, 
which is sensitive to the interference of unscattered and scattered waves, occurs over a region of spatial extent 
l/g'^T, which is the same as the mean free path for additional scatterings of the quark. 




FIG. 4. The generic diagram which will contribute to bremsstrahlung and inelastic pair annihilation 
at leading order. The solid outer lines represent charged particles whose momenta are hard but coUinear 
with the (hard, lightlike) incoming photon; These propagators can all be approximately on-shell simultane- 
ously. Therefore, self-energy resummation, including the imaginary part of the self-energy, is required. The 
cross-rungs represent soft gauge boson exchange with momenta ^ g^T^ which are restricted to satisfy 
Q ■ K ^ g^T^ in order to maintain the coUinearity condition. 



momentum parametrically large compared to gT. We will argue that we can also neglect 
any scattering process which changes the direction of the outgoing quark by an angle small 
compared to g. In particular, scattering off the g^T ultrasoft gauge boson background will 
not be important. This is because the 0{g^T) momentum exchanged in such a scattering 
is too small to affect the kinematics of the photon emission process. It is true that the 
scattered and unscattered waves also pick up an 0(1) color rotation (or a phase in an Abelian 
theory) from moving through the ultrasoft field over a distance 1/g'^T, but because they are 
overlapping and collinear, the color rotations they receive are the same to within 0{g), 
and so can not produce an 0(1) phase difference which would affect the photon emission 
amplitude. We will see how this physical argument arises in a diagrammatic analysis in 
subsection |1V E| . 

In this paper we show how to sum all processes like Fig. ^but with an arbitrary number of 
additional scatterings, working at all times to leading order in g. We show that the important 
exchange momenta are parametrically 0{gT); both harder and softer exchange momenta are 
irrelevant at leading order. Because the spatial extent of a gT scattering event is 0{l/gT), 
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while the mean free path is 0{l/g'^T), the scatterings can be considered sequential. Writing 
the photon emission rate in terms of a current-current correlator, or equivalently a photon 
self-energy, we show that the dominant set of diagrams for bremsstrahlung and inelastic pair 
annihilation are precisely ladder diagrams of the form shown in Fig. |. 



C. Results 



Let r^^^ denote the contribution to the photon emission rate per unit volume. 



dVdt 



;i-2) 



from bremsstrahlung and inelastic pair annihilation processes. The contribution Tl;^^^ should 
be added to the standard results p3| , p^ for 2 ^ 2 processes (Fig. 2) to obtain the total 



rate at (9(aEM«s)- We assume that the temperature is large compared to zero tempera- 
ture masses (or Aqcd), and neglect relative corrections suppressed by powers of g(T) or 

/ g^T"^. We also assume that the chemical potential for fermions is at most of order of the 
temperature, ^ < 0{T), but if scalar charge carriers are present, their chemical potential 
must be parametrically small compared to gT (since otherwise we would have to deal with 
Bose condensation of the scalars — something we wish to ignore). Finally, we require that 
the energy of the photon be parametrically large compared to g^T In g~^, so that the photon 
wavelength is much smaller than the large angle mean free path. 

Our result for the (unpolarized) differential emission rate dT]^^^ from a single complex 
scalar or Dirac fermion species in an equilibrium plasma at temperature T and chemical 
potential /i is given by 



d^k 



dp Qs Qe 



11^ /g^A(p||,fc)Re{2p^-f(p^;P||,A:)}. (1.3) 



d^p 



Here dp is the dimension of the quark's representation [A^^c in SU(Ac) gauge theory, or 3 for 
QCD], qs is the Abehan charge of the quark (2/3 for up type and —1/3 for down type), and 
k = |k|. The kinematic factor A{p\^,k) depends on the spin and statistics of the emitting 
particle. 



A{p\\,k) 



Mk+P\\) [l+ribipw)] 

2p\\ {p\\+k) 
nf{k+p\\) [1-nji 

2[p|| {p\i+k)Y 



+ iP\\+kf 



scalars; 



fermions, 



:i.4) 



with 



nb{p) 



exp[(3{p-fi)] 



nf{p) 



1 



exp[(3{p-fi)] + 1 



:i.5) 



the standard Bose and Fermi population functions, respectively, with fi the chemical poten- 
tial per quark (not per baryon). 

The quantity f{p±;p\\,k) appearing in the integrand of Eq. ( |1.3| ) is the solution to the 
following linear integral equation which sums the effects of multiple scatterings occurring 
during the photo-emission process. 
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2p_L = i6E i{p±;pu,k) + 



IT 



2 2/" d'^1± dq\\ dq^ ^ q 



X 



X 



T 

|q| 



[1 - (97lq|) 



212 



f(p±;pih^) - f(q+p±;Pii,^) 



:i.6) 



with Cp the quadratic Casimir for the quark [Cy = {N^—l)/2Nc = 4/3 in QCD], m-D the 
leading-order Debye mass, and 6E the difference in quasi-particle energies (defined as real 
parts of pole positions) between the state where the photon has been emitted and the state 
where it has not, 



6E = k^ + Ep sign(pii) - Ep+k sign(p||+fc) 



(1.7) 



For an SU(A^) gauge theory with A^^s complex scalars and A^f Dirac fermions in the funda- 
mental representation, the Debye mass is given by |31[ 



rriT 



6 



27r2 



(1.8) 



(where we have allowed a chemical potential fi for fermionic number but have taken the 
chemical potential for scalars to be negligible). In the denominators of Eq. ( [L.6| ), UtIQ) 
and Ili^{Q) are the standard transverse and longitudinal gauge boson self energies [given 
explicitly in Eq. ( |5.18| )]. 

In the final integral ( |1.3|) for the emission rate, the region p\\ > represents bremsstrah- 
lung emission by a particle of energy ~ A; + . In this region, the product of statistical 
distribution functions n{k+p\\)[l±n{p\\)] contained in A{p\\,k) obviously represent the pop- 
ulation function of an incoming quark of energy k+p\\ times the appropriate Bose enhance- 
ment or Pauli blocking factor for the outgoing quark of energy p\\. The region — |k| < p\\ < 
represents inelastic pair annihilation. Using the relations 



M-p) = -[1 + Mp)] > ^fi-p) = [1 - ^/(p)] 



:i.9) 



with n(p) = l/le'^^P^^^ =p 1] the appropriate anti-particle distribution function, the factor 
A{p\\, k) in this interval may be rewritten in the form 



A{p\\,k) = 



' ^&(^-biil)^fc(biil) 

2hl(A:-hl) ' 

^/(^-biil)%(biil) 

2[hl(A:-b|||)]2 



scalars; 



+ (fc— |j9|||)^ , fermions. 



:i.io) 



which displays the expected population functions for a particle and anti-particle to annihi- 
late, producing a photon of energy k. Finally, the region py < — |k| represents bremsstrah- 
lung emission by an anti-particle of energy | , and gives a contribution which is the same 
as the p|| > region except for switching n^j — > rihj- 

The integral equation ( |1.6| ) is similar in form to a linearized Boltzmann equation, but 
with the convective derivative replaced by ibE which represents a net phase accumulated due 
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to the energy difference between outgoing and incoming states. The integral in Eq. (|1.6|) can 
be interpreted as a hnearized colhsion integral, with the f (q+p^; py, fc) piece representing 
the gain term describing the scattering of particles into momentum p, and the f(p_|_;p||, k) 
piece the corresponding loss term describing the scattering of particles out of the mode with 
momentum p. An integral equation of this form was first derived by Migdal p9| , p0| in the 
context of energy loss of a fast particle traversing ordinary matter. 

The explicit form of the energy Ep of a hard quark with momentum |p| is given by 



p2 + 



IpI + 



iPl 



hi + 



2hl 



:i.ii) 



where the asymptotic thermal "mass'' 



:i.l2) 



(For scalar quarks with non-negligible quartic coupling, an additional term ~ AT^ also 
contributes.) As noted earlier, this asymptotic mass is the same for a fermion or a scalar, 
and the above approximation to the momentum dependence of the quasi-particle energy is, 
for hard excitations, accurate to 0{g^T) which is sufficient for our leading order analysis. 
The final form of Ep in Eq. ( |1.11| ) uses an expansion in p^ ^ which is justified because 
it is the region p^ = 0{g'^T'^) which gives the leading (in g) contribution to the emission 
rate; this expansion was already used in deriving the integral equation (|1.6| ). Substituting 
the explicit form of Ep into the definition (|1.7| ) gives 



6E 



pi + 



(1.13) 



for dispersion free, [Jon- shell photon emission, so that /cq = |k|. For off-shell photon emission 
(which can be relevant for lepton pair production), one must add to this expression the 
off-shell frequency shift k^—k. All of our analysis, including the justification that multiple 
scattering effects must be summed to determine correctly the leading photon emission rate, 
remains applicable to off-shell photons whose four-momenta satisfy /c| ~ g^T. When 
g'^T <C /c| -C T, our treatment is correct but unnecessarily complicated; scatterings 
can be ignored, as either pair annihilation (for > k) or bremsstrahlung (for /c > /c°) is 
kinematically allowed. For lA;*^— A;| > T, our treatment breaks down. 

The purpose of this paper is to present a derivation of Eqs. ( |1.3D and (|1.6|), based on a 
careful diagrammatic analysis showing that the ladder diagrams of Fig. ^, and only these 
diagrams, contribute to the leading order photo-emission rate. We begin with a somewhat 
more technical overview of the problem in Sec. Then we review our conventions and 



the diagrammatic analysis techniques we will need in Sec. 111. The power counting analysis 



''As noted earlier, the photon dispersion relation is also modified by medium effects, but the correction, 
kg — k'^ ~ e^r^, is smaller than the quark dispersion correction by 0(aEM/cks) <C 1. But if we were interested 
in the emissivity of an electromagnetic plasma, then photon dispersion corrections would need to be retained 
when evaluating 5E. 
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for the hard photon emission rate, leading to the identification of the diagrams shown in 
Fig. ^, is presented in Sec. This section also shows that the photo-emission rate is not 
sensitive to ultrasoft g^T scale physics, but only because of a cancellation between multiple 
diagrams containing ultrasoft gauge boson exchanges. The leading order ladder diagrams 
are summed in Sec. |V]. These two sections are the crux of the paper. We discuss various 
consistency checks, such as confirming that our approximations give a transverse current- 
current correlator and a finite photo-emission rate, in Sec. |VI[ 

For simplicity of presentation, the bulk of our analysis will be carried out using scalar 
quarks instead of fermions, and will focus on the case of hard, nearly on-shell photons; that 
is, photon momenta k = 0{T) with \K'^\ ^ g'^T'^- In the final section |V11| , we generalize 
the treatment to the case of fermions, and discuss the applicability of our analysis to softer 
and/or more highly off-shell photons. For on-shell photons, we find that the above results 
for the leading order photon emission rate remain valid provided the photon momentum is 
large compared to the large angle scattering rate, |k| 3> g'^Tlng'^. 



II. DIAGRAMMATIC APPROACH: OVERVIEW 

The differential unpolarized photon emission rate (per unit volume), at leading order in 
e^, is given by the well-known relation 

where K is the null photon 4-momentum with 3-momentum k and positive energy = 
|k| = k, and Wfj^i^{K) is the Wightman electromagnetic current- current correlator, 

W,.{K) = J d'x e-'^^ {j,{x)UO)) ■ (2.2) 

[We use a metric with ( — 1-++) signature.] As always, (■ ■ ■) denotes an expectation value 
in whatever density matrix is of interest, which in our case is a thermal ensemble describing 
the equilibrium plasma. 

The photon polarization basis vectors {e^^-)(k)} may be chosen to be unit spatial vectors 
orthogonal to k. Of course, gauge invariance ensures that the correlator is transverse, 

Wf^u{K) = 0, which implies that one may replace the sum of projections onto photon 
polarizations by g'^'^. However, for our analysis it will actually be more convenient to retain 
explicit transverse polarizations as in ( p.l|) . 

Note that the product of currents in the expectation value ( |2.2| ) is not time ordered, 
but is in the order shown. In other words, we need a Wightman current-current correlation 
function. In thermal equilibrium, this is proportional to the spectral density, or equivalently 
to the imaginary part of the retarded correlator, 

W,AK) = 2 K(A:°) + 1] lmD^:\K) , (2.3) 

with 

D^^\K) = I I d'xe~^^-Q{x') (b,(x),j.(0)]) , (2.4) 



11 



(A) (B) 

FIG. 5. Diagrams contributing to the current-current correlation function. The shaded blob in (A) 
represents the amputated four-point correlation function of the electrically charged fields which generate the 
current. Tadpole diagrams (B) make no contribution to the correlator of interest and may be ignored. 

and Q{x) the usual unit step function. [Since the current is charge neutral, the Bose 
distribution function in Eq. ( |2.3| ) does not include any chemical potential.] Because the 
current is bilinear in the fundamental fields, this correlator is determined by the 

connected four-point correlation function of the underlying (electrically) charged fields, as 
shown in Fig. ^(A). For theories containing charged scalar fields, the tadpole contribution to 
the retarded (or time-ordered) correlator, shown in Fig. ^(B), has no imaginary part. Hence 
it does not contribute to the Wightman function Wf^i, and may be ignored. 

The following discussion will, for simplicity, focus on the case of a theory with a single 
scalar "quark" with electric charge e. The generalization to fermions appears in Sec. |V11 A. 



We will label momenta as shown in Fig. ^(A), so that Pi is the momentum of the line 
entering the first current insertion, and P2 is the momentum of the line emerging from the 
second insertion. Given this, the photo-emission rate per unit volume may be expressed as 

= S^M ? ''"^'^"^ iK+2PMK+2P,)^ G.MPu -K-P,, K+P,, -P,) 

OEM rf^k f d^Pl d^P^ 



TT 



^ |k| / (2^ (2^ ^^'^ ■ ^^'^ ^+^^' -^'^ ' ('-^^ 



where we define to be the projection of a 4-momentum P onto the plane orthogonal to 
both K and k. In this expression, G'1122 denotes the Fourier transform (defined so that all 
momenta are incoming) of the scalar field four-point function 

Gu22{xuX2; yu I/2) = ( T {^^(y^) 0(yi) } T {^^(xa) (P{x,) }) , (2.6) 

where T denotes time ordering and T anti-time ordering. For our immediate purposes it is 
sufficient to note that the lowest order one loop diagram (where the shaded box is absent 
and Pi = P2) represents the disconnected part of G'1122 which is just a product of Wightman 
correlators, 

Gif^2ixuX2;yuy2) = (0^1/2) 0(xi)) (0(i/i) 0^(x2)) , (2.7) 

and that any equilibrium Wightman two-point function (Fourier transformed) is n{k^)+l 
times the corresponding spectral density [cf. Eq. (^.3])]. 

Throughout our analysis, we consider diagrams in which self energy insertions on internal 
lines have been resummed {i.e., skeleton diagrams). For sufficiently high temperature, so 
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+ CROSSINGS 



+ CROSSINGS 



FIG. 6. Non-coUinear contributions to the current-current correlator, and their interpretation as products 
of (quasi-particle) scattering amphtudes. The first diagram gives the squares of the amphtudes for individual 
scattering processes, while the second diagram generates the appropriate interference terms. 

that as{T) is small, this means that the spectral density associated with a line carrying 
momentum P has a quasiparticle peak at = 0{g'^T'^), together with a smooth "off- 
shell" background. The width in energy of the quasiparticle peak is 0{g'^T\ng^^), and its 
integrated spectral weight is 0(1), while the spectral weight of the off-shell background is 
0{g'^) |]T6|,^. [These estimates assume that the components of P are 0{T). The next 
section contains a more detailed description.] 

To understand which kinematic regimes can contribute to the leading order hard, on-shell 
photon emission rate [i.e., = with |k| ^ T], it will be helpful to consider separately 
the following regions.^ 

1. Non-coUinear photon emission. By this, we mean that the perpendicular components 
of momenta pi and p2 (relative to the photon momentum k) are both comparable to 
the temperature, so Pi^± and P2,± are each 0(T). 

2. Near-collinear photon emission. One or both of the momenta pi and p2 are nearly 
coUinear with k. Specifically, Pi^± and/or P2,± are 0{gT) or smaller. 

The leading contribution from region (1) is comparatively easy to analyze. For the lowest- 
order one loop diagram where P2 = Pi, the condition that Pi x be 0{T) implies that the 
propagators carrying momenta Pi and Pi+K cannot both be nearly on-shell. The dominant 
contribution arises when one propagator, say the one with momentum Pi, gives an on-shell 
spectral weight, while the other propagator with momentum K+Pi gives an off-shell spectral 



^For completeness, there is also the transition region in which Pi,± and/or P2,± are 0{g'^T) for some power 
f intermediate between zero and one, < < 1. Contributions from this regime could potentially give rise 
to logarithms of the form \n[T/{gT)] = \ng^^. But this region never dominates the contributions of both 
the non-coUinear and near-collinear regions by a power of g. Hence, for the purpose of identifying which 
classes of diagrams can contribute to the leading-order ph oton emiss ion rate, this transition region need not 
be separately considered. Nevertheless, our final results (|l.3| ) and (|l.6| ), together with existing treatments 
of the non-coUinear region, do correctly include contributions from this intermediate region. 
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weight (or vice- versa). Such an off-shell spectral weight may be viewed as a cut through a 
self-energy bubble, as shown in Fig. ^ This may be recognized as describing the lowest-order 
Compton and annihilation contributions to the photon production rate. The appearance of 
a bubble means that one must also consider other two loop contributions. As shown in 
Fig. ^, these generate the interference diagrams required to describe the complete lowest- 
order Compton and annihilation processes. These processes yield an 0{as) contribution to 
the current-current correlator, since there is always either a loop correction or an 0{g^) 
off-shell spectral weight involved. [For fermions, it turns out that there is also a logarithmic 
infrared enhancement as the off-shell line becomes soft {i.e., parallel and perpendicular 
components of momentum on one line both small compared to T).] As these contributions 
have been computed in detail in previous literature pSj^ , we will not discuss this piece of 
the photon emission rate in any detail. The contributions from region (1) for higher loop 
diagrams are suppressed by additional factors of gs.f\ 

The collinear region (2) is more difficult to analyze. The regime where Pi (or P2) is 
nearly collinear with K is phase space suppressed; just requiring the momenta to be in 
the right range provides a suppression, since both components of Pi^± must be 0{gT). 
Moreover, since Eq. ( |2.5|) contains a factor of Pi,± ■ P2,±, there is at least another 0{g) 
suppression. However, the contour for the frequency integration [over {pi)^] is nearly pinched 
between the on-shell poles of the Pi and K+Pi propagators, which leads to a "pinching pole" 
enhancement of the frequency integral by a factor of the inverse distance between these poles. 
These poles are displaced by 6E ~ 0{g'^T), which implies [as discussed in detail in section 
^Vl , c.f. Eq. ( [4.8|) ] that the frequency integral generates a l/g"^ which compensates the O^g"^) 
phase space suppression. Because the precise separation of the on-shell poles of the two 
propagators is relevant, the inclusion of correct self-energies on the propagators is of 0(1) 
importance. For the lowest-order one loop diagram. Pi ^ ■ P2,± = Pi± = 0{g'^T^) in this 
region, which means the contribution from this collinear region is 0(a;s), or equally important 
as the contribution from region (1). Moreover, higher loop diagrams in this collinear region 
may contain additional "pinching pole" enhancements which can compensate for the explicit 
vertex factors of g^ associated with each additional loop. This implies that an infinite class of 
diagrams may contribute to the leading-order emission rate. It has been argued and we 
will show below, that all uncrossed "ladder" exchanges of gauge bosons carrying momentum 
0{gT) contribute at the same order as the simple one-loop diagram. The next two sections 
show that such diagrams, and only such diagrams, are important, while Sec. will show 
how they may be summed to obtain Eqs. ( |1.3| ) and (|l.(j|). 



III. REAL-TIME THERMAL FIELD THEORY 

The correlation functions we are interested in have the general form 



^The power counting analysis of section FV may be easily adapted to verify this. As will be discussed in 
detail below, the only way higher loop diagrams could not be suppressed is if each additional loop involves 
either soft or collinear enhancements which are sufficient to cancel an explicit factor of g^. This is precisely 
what happens in the near-coUinear region (2), which will be analyzed in detail, but not in the non-coUinear 
region (1). In fact, the first corrections to the leading non-collinear contributions involve the soft sector and 
are suppressed by one power of g. 
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IN M \ I N M \ 

(t{ n o,]t{x{o?^)^Itc{ n o,,x{o,,]), (3.1) 

\ j=M+l i=l I \ j=M+l i=l I 

where T denotes time ordering and T anti-time ordering, so that, reading from right to 
left, there are a series of operators at ascending times followed by a series at descending 
times. The scalar field four-point function 6*1122 in Eq- ( |2.6| ) is precisely of this form. As 
is well-known, correlation functions of this type may be represented by functional integrals 
in which fields are defined on a complex time contour which doubly traverses the real time 



axis forward and back, and then runs into the lower half plane a distance —i(3 |]12[. In the 
second form of Eq. (|3.1| ), operators labeled 1 are to be viewed as insertions on the future- 
directed part of the contour, operators labeled 2 are insertions on the past-directed part of 
the contour, and Tc denotes contour ordering. 

All the usual rules of diagrammatic perturbation theory immediately generalize to this 
"closed-time-path" formalism. Because the expectation value implicitly contains factors 
of exp{—iHt) on the forward-directed branch of the time-contour, and exp{iHt) on the 
past-directed branch, interaction vertices on the 1 and 2 branches have opposite sign. For 
example, a interaction produces insertions of 

-tXj (cp\{x) - <p\{x)) , (3.2) 

where / dl^x denotes a conventional integral over Minkowski space. We will not review the 
basics for this closed-time-path formalism further, since it is well treated in several textbooks 



As noted in Refs. for calculational purposes it turns out to be quite helpful to 

make a change of basis from the upper and lower contour fields ^2} to the so-called r, a 
basis in which 

0r = ^ ^"^ , 0a = 01 - 02 , (3.3) 

and likewise for other fields.]^ The reasons for this will be noted shortly. Propagators in 
this basis are related to conventional retarded or advanced propagators, or to the spectral 
density, as follows.^ 

G„{P) = t [d^x e-^^^ {Tc {M^) 0j(O))) = z [I + n,ip')] p{P) , (3.4) 



-'^'^Interactions also generate insertions on the "stub" of the time contour extending into the lower half plane, 
which produce perturbative corrections to the density matrix describing the equilibrium thermal ensemble. 
These "Euclidean" contributions are not relevant to the current discussion. 

This basis was originally proposed by Keldysh [O, and i s sometim es referred to as the Keldysh basis. 



Properties of the r, a basis are discussed for instance m Refs. 136 3^j3j]. The notation is not uniform in the 



literature; another common notation is to write our (1,2) as (+, — ) and our (r, o) as (1,2) 

"'^Note that our propagators in the r, a basis are defined without any factors of 2, which differs from the 
convention of Also, note that the first r or a index on the propagator refers to the field 0, and the 

second index to 
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GUP) =ij d^x e-^^^ (Tc (M^) 01(0))) = G'Ret(P) , (3.5) 

GariP) =ljd''x e-*^^ {Tc (0a(x) 0t(O))) = GAdv(P) = [GRet(P)]* , (3.6) 

Gaa{P) ^ijd'x e-'^^ {Tc {Ux) 01(0)) ) = . (3.7) 

Here and throughout we have suppressed writing group indices on propagators. The retarded 
and advanced propagators are0 



GRet(P) =z j d^xe-^^^e(x°)([0(x),0t(O)]), (3.8) 
GAdv(P) = ~ij d^x e-*^^ 0(-x°) ([0(x), 0t(o)]) , (3.9) 
while the spectral density 

p{P)^ /ci^xe-^^-([0(x),0t(o)]) 

= I {GAdAP) - GkAP)) = 21mGnet{P) ■ (3.10) 

The spectral density p{P) is real, odd in p^, even in p, and positive for positive Hence 
Grr{P) is symmetric in both p and p^, and is pure imaginary with positive imaginary 
part. Up to a factor of i, G„- is the same as the ordering- averaged two point function, 
G„iP) = tJd^xe-'P^ (i{0(x),0t(o)}). 

The statistical factor | + ni,ijp) = |coth(p°/2T) appearing in relation (|3.4| ) applies to 
Bose fields only; for fermionic fields it is replaced by | — nf{p^) = ^ tanh(p°/2T). Note 
that both I + nfc(p°) and | — rifijp) are odd functions of Also, for fermionic fields the 
commutators in Eqs. (p.8|), ( |0| ) and ( p.lO| ) are replaced by anticommutators. The symmetry 
and reality properties of the spectral density are the same for Bose and Fermi fields. 

Because interaction vertices are necessarily odd under the interchange of 0i and 02, as 
illustrated by Eq. (|3.2|) , in the r, a basis all vertices must involve an odd number of a indices. 
In particular, there are no pure r vertices. This fact, together with the vanishing of Gaa, 
helps make the r, a basis particularly convenient. In the following, we will need the explicit 
form of the vertex characterizing the scalar-field electromagnetic interaction. If the scalar 
lines leading into and out of the vertex carry momenta P and P', respectively, then the 
interaction with a single gauge field is 
P\ 

W = ig [P+P']^ {A't $1 $1 - $^ $2} 



P 



[P+P\ {A^ ($t + $t $^) + ($1 + 1 $1 $,) } , (3.11) 



where $ carries momentum P, $^ carries —P', and carries P'—P. 

In the free field limit, a retarded scalar propagator has the standard form 



^^Because is a complex field, GRct{—P) ^ GAdv{P), unlike the case of a real field; rather GRcti—P) is 
the advanced propagator which differs from Eq. ( |3.9|) by interchanging (p and (p^ . This implies that changing 
the sign of the four- momentum does not convert TT^a iuto Gar- 
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G^-t(P) = [p2 + - + itf] ' = (G^''^{P)X . (3.12) 

free L J \ / free 

In the interacting theory, a non-trivial self-energy S(P) will also appear in the denominator. 
At finite temperature, the self-energy has a non-zero imaginary part for all nonzero real 
This retarded self-energy is obtained from the Euclidean self-energy by continuing to real 
frequencies from positive imaginary Matsubara frequencies. 

For hard momentum and energy, Ip],/)*^ = 0{T), the retarded propagator G^^^{P) for any 
basic field (scalars, fermions, or gauge bosons) has a self-energy whose real and imaginary 
parts are both Olg^T"^)!^ This applies both on- and off-shell.0 Therefore, the spectral 
weight p{P) is 0{g'^/T'^) if the momentum P is far off-shell [P^ = ©(T^)]. The propagating 
pole of the retarded propagator acquires an 0{g'^T) imaginary part in frequency space, so for 
near on-shell hard momenta [P = 0(T) with P^ = 0{g^T'^)] the spectral density resembles 
a Lorentzian whose peak value is 0(1/ g^T'^) and whose width is 0{g'^T). And for hard 
energies, Grr{P) has the same characteristic size as the spectral density p{P) (either on- or 
off-shell) since the Bose (or Fermi) distribution functions are 0(1). 

For soft momentum and energy, |p|,p° = 0{gT), the real part of the retarded self- 
energy for scalars, gauge bosons, and fermions remains 0{g^T^).f^ The imaginary part is 
0{g^T^), except for gauge bosons and fermions with spacelike momentum P, in which case 
Landau damping produces an 0{g^T'^) imaginary self energy PHig. Note that for soft 



excitations with energies of order gT, an 0{g^T'^) on-shell imaginary self-energy implies 
that the quasiparticle width is 0{g'^T) (up to logs), just as it is for hard excitations. 

Consequently, for generic soft momenta (and energy), the retarded propagator for any 
excitation is 0{1/ g'^T'^). For spacelike gauge bosons, the corresponding spectral density is 
0(1/ g'^T'^), while Grr{P) is Oll/g^T"^) due to the T/p° ~ 1/g enhancement from the Bose 
distribution function. For soft scalars or timelike gauge bosons, the additional factor of g 
in the 0{g^T'^) imaginary self-energy implies that the spectral density is only 0(1/ gT"^) so 
Grr{P) is 0{1/ g'^T'^). These estimates hold provided the energy is off-shell (relative to the 
pole position defined by the real part of the self-energy) by 0{gT). If the energy is within 
0{g^T) of being on-shell, then the nearly on-shell retarded propagator and the spectral 
density are both 0{l/g^T'^), and Grr for scalars or gauge bosons is 0{l/g'^T'^). The above 
estimates are summarized in Table |. 



-"^^In the case of (massless) fermions, whose propagators have the form ^ = ( ^)/(P+E)^, 

our statements about the size of the "propagator" are statements about 1/(P+E)^, and statements about 
the "self-energy" are really estimates for the deviation of the rationalized denominator from its free value; 
that is, for 2P- S + E^. 

-'^^The imaginary part of the on-shcU self-energy has a logarithmic enhancement M] and is 0{ g^T^ Ing^^) 
due to soft scattering with momentum transfers extending from 0{gT) down to 0{g^), see Eq. ( |5.15 ). Since 



the goal our power-counting analysis is to classify contributions according to powers of such logarithmic 
effec ts ar e not relevant and will be ignored, until we reach the discussion of ultrasoft contributions 



Sec. [VE 



i^These estimates for soft momenta assume vanishing or negligible chemical potential. It will turn out 
that we actually only need the results for soft spacelike gauge bosons (which are unaffected by the chemical 
potential); soft charged fields will not play a role in our leading-order analysis. 
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Hard momentum 


Soft momentum 


Quantity 


off-shell 


on-shell 


spacelike 
gauge boson 


on-shell 
boson 


other 
boson 


G'Rct(-P) 






i/{g'T') 


1/ig'T') 


Uig'T') 


P{P) 






yig'T') 


i/ig'T') 


UigT') 


Grr{P) 




yig'T') 


i/{g'T') 


i/ig'T') 


i/ig'T') 



TABLE I. Characteristic sizes of various propagators in different momentum regions. "Hard, off-shell" 
means that |p|,p° = 0(T) with = 0{T'^), "Hard, on-shell" means |p|,p° = 0{T) with ^ 0{g'^T'^), 
and "Soft" means IpIjP*^ — 0{gT). For hard momenta, the estimates shown apply to scalar, fermion, or 
gauge boson propagators. For soft momenta, "on-shell" means that the energy is within 0{g^T) of the 
appropriate scalar or gauge boson quasiparticle peak, and "other boson" denotes either a scalar propagator 
with arbitrary soft 4-momentum P, or a gauge boson with soft timelike P, provided P is outside the 
respective on-shell regions. For soft fermions, all results are the same as for soft gauge bosons except that 
Grr is smaller by a factor of g, as there is no 1/g enhancement from the Fermi distribution function. The 
estimates for soft scalars assume negligible chemical potential. 



IV. POWER COUNTING 




FIG. 7. The important class of diagrams, in which both photon lines connect to the same quark loop. 
The blob in the middle connects gluon lines in an arbitrary fashion. It need not be connected, and may 
contain arbitrarily many non-Abelian vertices, additional quark loops, and gluon propagators. 



A. Presence of pinching poles 

We want to evaluate the contribution to the photon emission rate ( p.5| ) arising from the 
kinematic region Pi^± = 0{gT). Diagrams which can contribute at leading order in the 
electromagnetic coupling e, but any order in Qs, include those of the general form shown in 
Fig. 1^ in which both photons connect to the same quark loop. Any number of gluon lines, 
which are themselves connected in an arbitrary fashion, may also attach to the quark loop. 
(Recall, however, that all propagators are to be regarded as containing full self-energies; 
only skeleton diagrams need be considered in this discussion, except for a portion of section 
[IV E| below.) In addition to these diagrams, there are also diagrams in which the two photon 
lines connect to different quark loops which are themselves connected by two or more gluon 
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FIG. 8. Simplest non-vanishing diagram in which the photon Hnes connect to different quarlc loops, 
which are themselves connected by two or more gluon lines. This class of diagrams will be seen to be 
sub-leading compared to the leading diagrams where both photons connect to the same quark loop. If the 
quark chemical potential vanishes, then charge conjugation invariance implies that there must be three or 
more giuons connecting the two quark loops. 




FIG. 9. Single gluon exchange diagram. 

lines, as illustrated in Fig. |[|^ The bulk of our discussion will focus on the first category of 
diagrams, but we will return to this second category in the next subsection. 

It will turn out that a restricted but infinite class of higher loop diagrams contribute at 
the same order as the lowest order one-loop graph. One may understand why as follows. 
First, consider a single gluon exchange, as shown in Fig. As noted earlier [cf. Fig. ^, 
this diagram generates the interference part of the leading 0{as) non-collinear contribution. 
But focus instead on the near-collinear region where P± is 0{gT) and the exchanged gluon 
momentum Q is soft, so that g° ~ gT and |q| ~ gT (and P+Q is also nearly collinear). We 
wish to compare the contribution from this near-collinear region to the 0{as) contribution 
from the non-collinear region. 

For exactly the same reasons discussed at the end of Section ||, the dp^ frequency in- 
tegration is dominated by the nearly pinching poles of the P and K+P propagators which 
happens when P is almost null, P^ = 0{g'^T'^), and nearly collinear with K, so that K+P 
is also almost null, {K+Py = 0{g'^T'^). Exactly the same reasoning implies that the domi- 
nant contribution to the dq^ frequency integral will occur when the poles of the P+Q and 
K+P+Q propagators nearly pinch, which happens when (P+Q)"^ and (K+P+Q)"^ are both 
0{g^T'^). These conditions are equivalent to the statement that Q^, Q ■ K, and Q ■ P are 



single gluon connecting the two quark loops is forbidden by (non-Abelian) gauge invariance. If the 
chemical potential of the quark is negligible, then charge conjugation invariance implies that there must be 
three or more giuons connecting the two quark loops; and if in addition the gauge group is real or pseudo-real 
(for instance, SU(2)), then a loop with one photon and any number of giuons vanishes. 
At higher orders in Oem, diagrams with multiple quark loops separated by one or more photon lines are 
also possible. One-particle reducible diagrams with single photons connecting quark loops represent re- 
absorption corrections; they can be easily included by interpreting our result as a computation of emissivity 
rather than the emission rate. For simplicity, we neglect all higher order in ckem corrections. 
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all 0{g^T'^). Now, the gauge boson propagator G^"{Q) has its indices contracted against 
the two 0{T) momenta {2K+2P+Q) and (2P+Q), which we have just argued are approxi- 
mately lightlike and collinear. If G^^ were the free propagator, it would behave like g^'^ and 
hence produce a g"^ suppression from contracting nearly collinear and lightlike 4-momenta. 
However, for a soft exchange momentum Q ~ gT, the thermal self-energy in the gauge 
boson propagator represents an 0(1) modification. Moreover, the transverse and longitudi- 
nal self-energies differ significantly, which means that the time and space components of the 
4-momenta which G^^ contracts contribute differently. Hence, there is no such near cancella- 
tion. As noted in the previous section, ~ l/g^T"^ for soft spacelike momenta. Therefore, 
the gluon propagator contracted between {2K+2P+Q) and {2P+Q) can be 0{l/g^). 

In the near-collinear region, there is a g^ phase space suppression from the soft (i^q inte- 
gration, and a g^ phase space suppression from the requirement that Pi be 0{gT) (so that 
p is nearly collinear with k). Finally there is the explicit g^ from gluon exchange vertices, 
and a g^ from the Pi ■ {P+Q)± factor associated with the photon vertices. Combining the 
pieces, one has a. g^ x g^ phase space suppression, g^ from vertex factors, two pinching-pole 
enhancement factors contributing (7^^ x g^"^, and a g~^ soft gluon propagator. The net result 
is an 0{g^) contribution, exactly the same as the leading-order non-collinear contribution. 

Moreover, since the momentum P+Q satisfies the same coUinearity conditions as P, 
adding additional soft gauge boson "cross-rungs" will also give relative 0(1) corrections as 
each additional rung adds an explicit g^ x g^ suppression from the new vertices plus soft phase 
space for the new exchange momentum, which is just balanced by a g^^ x g~^ enhancement 
from one more pinching-pole frequency integral plus another soft gluon propagator. This 
will be shown in greater detail in the next section. Hence, at a minimum, all ladder diagrams 
will contribute to the leading order emission rate. 

Physically, these gauge boson exchange corrections are important because the lowest 
order one loop diagram is dominated by a narrow range of frequencies which is 0{g^T) in 
width. That means that the time scale for the emission process is ~ 1/g'^T, which is the 
same (up to logarithms) as the mean free path for small- angle scattering. Hence, 0(1) effects 
from such scatterings are to be expected. 

A similar need to resum ladder diagrams is also encountered when computing transport 
coefficients in many other contexts including the textbook case of conductivity in ordinary 
metals (see, for example, Ref. |^8|] and Sec. 39.2 of Ref. p^), viscosity in relativistic scalar 
field theory or low frequency electromagnetic response in QCD ^T|. In all these appli- 



cations, the presence of nearly pinching pairs of poles also necessitates the summation of an 
infinite set of ladder-like diagrams. However, because we are considering hard, lightlike ex- 
ternal momentum, rather than very soft, the appropriate power counting analysis differs in 
detail, and cannot simply be taken over from previous work. Although the relevant diagram 
topologies which must be summed are superficially the same, the kinematic regions of the 
diagrams which are important at leading order prove to be rather different. In particular, 
we will find that only soft momentum flowing through the cross-rungs will contribute in our 
application of photo-emission; this is not true for transport coefficients in either scalar field 
theory or gauge theories. 
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FIG. 10. A lowest-order example of a diagram containing a three point non-Abelian vertex. This, and 
all other diagrams containing three-point vertices, are subleading. 

B. Absence of non-Abelian vertices and multiple quark loops 

For "nice" problems in thermal field theory, in which T and gT are the only relevant 
momentum or frequency scales, perturbation theory using HTL resummed propagators and 
vertices yields an asymptotic expansion in powers of g (not g"^), with each order resulting 
from only a finite set of diagrams. This simple diagrammatic structure breaks down whenever 
dependence on time scales comparable to the small-angle mean free scattering time, which 
scales as l/{g'^T) (up to logs), becomes relevant. It also breaks down whenever dependence 
on gauge field fluctuations with spatial scales of l/{g'^T) or larger becomes relevant, because 
the statistics (and dynamics) of such fluctuations is non-perturbative. 

The first condition, dependence on {g'^T)^^ time scales, is the reason why higher loop 
diagrams can contribute to the leading-order hard photon emission rate. As noted earlier, 
this dependence arises from the nearly pinching poles in frequency integrals, which reflect the 
presence of virtual intermediate states that propagate for times of order {g'^T)~^. However, 
the only higher loop diagrams which can contribute to the leading-order emission rate are 
those in which each additional loop brings with it one more pinching-pole enhancement 
factor. In particular, diagrams in which the "blob" in Fig. |^ contains vertices and loops 
are suppressed. (Recall that self-energy insertions on all lines are implicitly resummed. We 
are not saying that self-energy insertions are sub-dominant on either the quark propagators 
forming the perimeter of the loop, or on the soft gluon propagators.) 

As an example, consider the diagram shown in Fig. and compare its contribution to 
the one-gluon exchange diagram of Fig. |^. Consider the case where all gauge boson lines are 
soft [0{gT) momenta]. There is an additional g"^ suppression from the explicit couplings on 
two more gluon vertices and a further g^ suppression from phase space, from one more soft 
spatial momentum integral. Since the three point vertex is a derivative interaction and all 
the lines entering it are soft, the vertex contributes an extra 0{g) suppression. That makes 
a g^ suppression so far. Two of the gauge boson propagators can be rr propagators, but 
the three point vertex must have one a index, so one propagator is not. For soft spacelike 
momenta, the rr propagators are each 0(1/ g^T^), while the soft ra (or ar) propagator is 
0(1/ g'^T^); hence the two additional gauge boson propagators make at most a g~^ contri- 
bution. There is one new frequency integral and one new scalar propagator. If the new 
scalar propagator is on-shell it will contribute l/g"^, but this requires an 0{g'^) restriction on 
the new frequency integration, for a net g^ contribution. Alternately we could say that one 
frequency integration involves a single off-shell scalar propagator and again contributes g^. 
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Either way, totaling it up, one finds that the diagram is 0{g) suppressed.[^ Had the gluon 
propagators met at a loop rather than a vertex, the dominant contribution would be from 
the HTL part of the loop, which is of the same order as the tree vertex, so the counting 
would be unchanged. 

The diagram in Fig. |10| is even more suppressed when a hard momentum flows through 
one or more of the gluons. Then, except in a phase space restricted range of momentum, 
there is not a second pair of pinching poles, and the diagram is suppressed, just as the 
most naive counting would suggest. The contribution of hard but collinear gluons is also 
suppressed (though by less than g"^) for the same reasons discussed in the next subsection. 
The case where the large loop momentum flows through the upper two gluons, and the upper 
middle scalar propagator is soft, is also suppressed, as there is no new pinching pole (see 
Sec. [IV D| ), and the gluon propagators C^" contract nearly null, collinear momenta, leading 
to an additional suppression. 

More or higher order vertices lead to further suppression, because they fail to introduce 
new pinching poles. For instance, adding another gauge boson line to Fig. |10| and replacing 
the three point vertex by a four point vertex brings in one more g^ phase space suppression 
(in the soft momentum region), compensated by a l/^f^ from the additional soft gauge boson 
propagator. The explicit g"^ on the four point vertex counts the same as the g of the three 
point vertex together with the extra g from the derivative coupling in the three point case. 
There is one more explicit g from the new scalar-gauge vertex, a new frequency integral and 
one more scalar propagator, but no additional pinching pole pairs. Hence, relative to the 
three point vertex, the four point vertex is suppressed by one more power of g. Essentially 
the same arguments rule out diagrams containing a "seagull" scalar-gauge field vertex, such 
as the one illustrated in Fig. ^, or diagrams where one of the photons emerges from a 4 point 
photo-gluon vertex, also illustrated by Fig. |^. 

Diagrams of the form illustrated in Fig. ^ are also suppressed. If the gluons are hard and 
collinear, and the quark "crossbars" (the vertical lines in the figure) are soft, the diagram is 
suppressed by g^ with respect to the leading order processes; a factor of g^ arises because the 
size of Grr in the soft, space-like region is g smaller for a quark (either scalar or fermionic) 
than for a gluon, and a factor of g^ arises from each gluon because they contract (almost) 
lightlike collinear momenta (see next subsection). The case where the quark crossbars are 
hard is analogous to the cases discussed in the next subsection. Also note that, in a C 
invariant theory at zero chemical potential, this diagram vanishes; three or more gluons 
must run between quark loops to provide a nonzero result, and such diagrams are even more 
suppressed. 

Hence, neglecting for the moment the possibility of "ultrasoft" gauge boson lines carry- 
ing momenta of order g'^T, the only diagrams we need consider are those where all gluon 
propagators attach at both ends to the "main" quark loop. That is, we need now only 
consider diagrams like that shown in Fig. |11|. We will return to ultrasoft momentum bosons 
later, as they require a more involved analysis. 



^^As noted in Table |, the sizes of soft gauge boson propagators are even larger, by one power of 1/g, for 
on-shell momenta as compared to soft spacelike momenta. But in the frequency region where a gauge boson 
propagator becomes on-shell, a scalar propagator leaves the mass shell, so there is no actual enhancement. 
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FIG. 11. An example of the remaining class of diagrams which must be considered, in which all gauge 
boson lines connect (at both ends) to the main quark loop. 

C. Momenta harder than gT 

The power counting of hard gluon exchange is relatively simple, so we will address it 
next. Consider Fig. ^ once again, but now when the momentum Q is hard, Q ^ T. For 
generic momentum q, the mass shells of Q, K+P+Q, and P+Q do not occur at the same 
value of so the g° frequency integral over the three propagators is 0(1). The explicit 
from the added vertices therefore suppresses the contribution. Consider instead the case of 
collinear momentum q, that is, q^ ~ gT. In this case, all scalar propagators and the gauge 
boson propagator can go on-shell simultaneously. This kinematic region is g"^ suppressed, 
as both components of q_L must be 0{gT). It also has the explicit g"^ from the vertices. At 
values of the frequency where is far off shell, there is no enhancement, so this region is 
suppressed. The frequency range of g° which brings all propagators on-shell is 0{g'^T) wide, 
and there are three new l/g"^ propagators, giving 1/g'^. But in this regime, the gauge boson 
propagator is contracting near-collinear momenta. For hard momenta, the Lorentz structure 
of G^"" is approximately gf"" (in Feynman gauge), but g^"'{K+2P+Q)^{2P+Q)„ ~ g'^T'^. 
Therefore, this hard but collinear region is also g^ suppressed.P^ 

Physically, the suppression of 0{T) scattering momenta comes about because the time 
scale involved in the emission process is 1/g'^T, but the mean free path between large angle 
scatterings is 0{1/ g'^T) (up to logarithms). Hence, the chance of a hard scattering occurring 
during the photon emission process is 0{g^). 

D. Crossed rungs and soft vertex corrections 

The next task is to rule out crossed rungs. Unlike the case of scalar field theory , this 
is not just a matter of counting the parametric size of each momentum integral. Consider 
the diagram drawn in Fig. It is not of the form of the diagrams in Fig. ^, and we claim it 
does not contribute at leading order. Yet, when all gluon momenta are 0{gT), simple power 
counting arguments suggest the diagram should contribute. Each gauge boson propagator 



-'^^The assumed collinearity conditions also imply that Q- {2P+Q) and Q- {K+2P+Q) are 0{g^T^). Hence, 
the conclusion that the hard but collinear region is 0{g'^) suppressed does not rely on the particular choice 
of Feynman gauge. 
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introduces a from its vertices and a of phase space suppression. But each gauge boson 
propagator, for soft spacehke momenta, has a spectral weight times thermal distribution 
function which is 0{1/ g^T"^). Further, each line has an 0{gT) transverse component, so 
they do not disturb the collinearity conditions on the hard momenta carried by the quarks. 
In a restricted frequency range, in which each loop frequency is tuned within an accuracy of 
g^T, the quark lines are on-shell (within g'^T in energy) simultaneously, providing al/g^ per 
propagator. Naively adding powers of g, this diagram should contribute at leading order in 
g^. Nevertheless, we will see that when the frequency integrations are actually carried out, 
the diagram proves to be suppressed. 

The analysis in this subsection is applicable when the momenta carried by gluons are all 
much larger than the g'^T ultrasoft scale; these arguments will not be applicable to ultrasoft 
gluons. However, ultrasoft gluons carry so little momentum that if any diagram with purely 
gT scale gluon lines is "clothed" by adding additional ultrasoft gluons, this will not affect any 
of the arguments of this subsection concerning gT scale gluons. Hence it is safe to postpone 
treating ultrasoft gluons, and focus first on the analysis of gT scale gluon exchanges. 

Crossed diagrams are suppressed, but uncrossed ladders are not, because of the details 
of the frequency integrations. The analysis turns out to be simplest in the r, a basis. But 
to use this basis, we must first transform 6*1122 into a linear combination of r, a basis four- 
point functions. By inverting the relation between r, a and 1, 2 bases given in Eq. (|3.3| ), one 
may easily express G1122 as a linear combination of the 16 possible four point r, a correlation 
functions. However, not all of the 4-point functions are independent in the r, a basis, and the 
relations between them have been worked out recently by Wang and Heinz [^. They show 



that (71122 can be expressed as a linear combination of just 7 of the 16 r, a basis four-point 
functions and their charge conjugates,^ 

6*1122 = Oil Gaarr + Ct2 Gaaar + «3 Gaara + Ct4 Garaa + «5 Graaa + OiQ Garra + «7 Garar 

+ Pi G^aj.j. + (32 G^^^^ + 13^ Gaara + l^i ^araa + G^^^^ + (3^ G^^^^ + (3-j G ^^^^ , (4.1) 

where all correlation functions have arguments ("Pi, ^2, ^3, "^4)- Here, 

Ga,...A4(Pi,P2,P3,P4) = / nci'a;,e-^^-^'(rc«(a;4)0A3(x3)0L(x2)0A,(xi)}), (4.2) 

■' i=i 

while G\^...x^ denotes the charge conjugated correlator, 

Ga,..A4 (Pi, 7^2,^3,^4) = J Ud'^i e-'''""'(rc{K(x4)0i3(x3)0A.(x2)0L(xi)} 



Gx,x,x,x,iV2,VuV^,V3). (4.3) 



^"Actually, Wang and Heinz 
to complex fields. Using only 



only discuss the case of a real scalar field, while we need the generalization 
PT invariance plus the KMS conditions satisfied by thermal equilibrium 
correlation functions, one may show that the appropriate generalization of the results of Ref. |^^, for bosonic 
fields, merely replaces the complex conjugate G* ^...^^ appearing in Ref. |34| by the complex conjugate of 
the charge conjugated correlator. If the equilibrium ensemble is charge conjugation invariant (which implies 
vanishing chemical potentials) then there is no difference. 
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If this index / is 'a 




This one , must be 'r' 



FIG. 12. Assignment of r, a indices in a ladder diagram. Because each vertex must have an a index, the 
gauge rungs must be rr, and the aa propagator vanishes, an a index at the beginning of a ladder "side-rail" 
forces alternating r, a assignments for the rest of the rail, with the other end necessarily ending in r. 

If the chemical potentials vanish, so the plasma is charge conjugation invariant, then 
Ga,...a4 (^1,^2,^3,^4) = Gx,...x,{Vi,V2,V^,Va). With our conventions^ 



ai = n[Vi] n[V2] 



Pi 



(4.4) 



with n[V] = nh{V^) shorthand for the indicated Bose distribution function. The remain- 
ing coefficients 02 ... ay and (32 - ■ ■ (3j involve similar combinations of distribution functions 
whose explicit form can be found in Ref. 



I , but will not be needed for our analysis as none 
of these other terms will turn out to contribute at leading order. 

The incoming momenta Pi • • • P4 are related to our previous choice of momentum labels 
[as shown in Fig. | through (Pi, P2, ^3, ^^4) = (Pi, -isT-Pi, i^+P2, -^2)- Since the mo- 
menta Pi and P2 are hard but all exchange momenta are soft, at leading order P4 ^ —Pi 
and P3 ^ — P2, in which case. 



/?i ^n[Pi]n[p2] = ai. 



(4.5) 



Consequently, Gaaw aarr 



and G^^j.^. enter with the same coefficient. The contribution of 
Gaarr the cmissiou rate (|2.5| ) is just the complex conjugate of the Gaarr contribution 
(even with non-zero chemical potentials), as may be seen by using the relation ( [4.3| ) to 
rewrite Gaarr{Pi, -K-Pi, K+P2, -P2) as Gaarr{-K-Pi, Pi, -P2, K+P2), making a shift 
Pi — s> —K—Pi, P2 —K—P2 in the integration momenta in Eq. ( p.5| ), and noticing that 
this shift has no effect on the other factors in the integrand. As a result, it is only the real 
part of GaaTr which contributes to the emission rate. 

As we have already seen, each new gluon propagator introduces a from vertices and a 
phase space suppression, as well as a Ijg"^ from the pinching-pole frequency integration 
involving the two new nearly on-shell scalar propagators. To contribute at leading order, 
the gauge boson propagator must be 0{g~^), which means it must be an rr propagator 
(see Table |). As shown by Eq. ( |3.11| ), every vertex must have one r and one a scalar 



^^Our conventions differ from jSj 
the correlator. 



their G is i T' 



times our G, where rir is the number of r indices on 
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index. Further, Gaa vanishes. According to Eq. (^J]), we need only consider diagrams 
for four-point functions with at least two a indices. Consider a propagator emerging from 
a photon vertex, with an a assignment at the vertex. The other end must be an r, as 
Gaa = 0. Since the gluonic vertex it goes into must have an a, and the gluon uses the 
other r, the next scalar propagator emerging from the vertex must begin with an a index. 
Regardless of the arrangement of the gauge propagators (crossed, uncrossed, etc.), the scalar 
line which enters the photon vertex on the other side has an a index on the gluonic vertex it 
leaves, and must therefore enter the photon vertex with an r index, as illustrated in Fig. 
Given the convention for ordering indices on the four-point function shown in Eq. ( |2.6| ), this 
demonstrates that only the Gaarr and Garar terms in Eq. (^?1| ) can contribute at leading 
order. 

Consider the two quark propagators leaving the first photon 



Return now to Fig. 11 



vertex, and regard P as the loop momentum circulating around the left-most loop of the 
diagram. Fix for the moment the three spatial components of P, and consider the integration 
over the frequency p*^ . As far as this frequency integral is concerned, the essential part of 
the integrand is the product of scalar propagators with momenta P and K+P; all other p° 
dependence is sufficiently weak that it may be neglected over the relevant frequency interval 
which is only 0{g^T) wide. This will be discussed in greater detail momentarily. If the 
photon vertex is an aa vertex, then one has 



J ^ Gar{K + P) Gra{P) = / ^ G^a^K + P) Gnet{P) 

1 n 



(4.6) 
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(p° + 1 rp)2 - E^] [(/ + fco - i r 



2 ^ p+k)^ - -£^p+k 



Here, Ep is the on-shell quasiparticle energy, Fp = lmT,{Ep,p)/ Ep is the quasiparticle 
decay width generated by the imaginary part of the on-shell self-energy, and the Lorentzian 
approximation to the propagators is valid for frequencies near the quasiparticle poles. [In 
thermal field theory rp/2 is often called the damping rate.] The integrand has four poles, 
two above and two below the real axis. The pole locations are 



p 



±Ep - iTp/2 , 



p 



-e ± Ep+k + ^rp+k/2 . 



(4.7) 



Since p and k are nearly coUinear, -Ep+k — Ep ± |k| (with the ± depending on whether 
p and k are aligned or anti-aligned). In either case, two pole positions almost coincide, at 



approximately p 



= p ■ k. One is from the retarded propagator and one is from the 



advanced propagator, so they are on opposite sides of the contour. That is, the poles nearly 
pinch the integration contour. The integral is dominated by the region near the pinching 
poles and one finds 



-^GaUK + P) GKotiP) 



1 



4=p\\{p\\+k) |(rp+rp+k) + «<5^ 



(4.8) 



where we have defined 



6E = 



Ep sign(p|| 



P\\ 



^p+kSign(p||+A;) - p\\ - 



(4.9) 
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as in Eq. ( |1.7| ), and approximated sign(p||) sign(p||+fc) EpEp^^ by P\\{p\\+k) in the denomi- 
nator of the result ([4.8| ) , since the difference is subleading in g in the relevant domain where 
p\\ and p\\+k are 0(T) while pj_ is 0{gT). Note that 6E and Fp are both 0{g^T) in this 
domain [assuming, as we do throughout this discussion, that the photon is off-shell by at 
most a comparable amount, k—k^ = 0{g'^T)]. Consequently, the width of the frequency 
interval which provides the dominant contribution to the integral is of order g^T. The fac- 
tor of ^(rp-t-Fp+k) + iSE in the denominator of ( [4 .81 ) is nothing other than {i times) the 
separation between the nearly pinching poles in the complex frequency plane. It should be 
borne in mind that 6E is a function of both p and the photon four-momentum K, even 
though this is not indicated explicitly. For hard momentum, the quasiparticle decay width 
Fp is essentially constant in p and differs from the asymptotic value F = limp_+oo Fp only 
by an amount which is subleading in g. As this will have no bearing on any of our results, 
henceforth we will generally simplify our expressions by omitting the spatial momentum 
dependence of Fp. 

The analysis leading to Eq. ( [4.8|) ignored the frequency dependence of the other parts 
of the integrand such as the gluon propagator, and also ignored frequency dependence in 
the imaginary part of the self-energy as one moves away from its on-shell value. But such 
frequency dependence has a characteristic scale of gT or larger, and so is negligible over the 
parametrically small 0{g^T) frequency region which dominates the integral. 

In contrast, if the photon vertex has one r and one a index, as occurs for Garar, then the 
integral we must consider is 

I^GnetiP + K)GnetiP), (4.10) 

and this has the same pole positions, except that the signs of the imaginary parts are all now 
the same. Therefore the contour for the frequency integration may be deformed away from 
these poles [by a distance large compared to the 0{g^T) displacement of the poles from the 
real axis], which implies that the integral has no l/g"^ enhancement, unlike the pinching-pole 
result n) . 

Therefore, of the various different r, a basis correlators appearing in Eq. {^A\) , only Gaarr 
(and its charge conjugate) contributes at leading order; Garar and all the others may be 
dropped. This substantial simplification is what makes the r, a representation convenient 
for this problem. A similar reduction has been noted recently for shear viscosity in scalar 
field theory H. 

Although it is not strictly necessary, it may be helpful at this point to verify explicitly 
that ladder diagrams with uncrossed rungs really do contribute at leading order — there is 
no unforeseen cancellation which would make their contribution subleading. Consider the 



diagram fragment shown in Fig. 13, and integrate over the frequency flowing around the 



^^This assertion is actually a bit too cavalier, because Fp is logarithmically sensitive to ultra-soft g^T scale 
physics and the g^T scale contributions to Fp do change significantly in an 0{g^T) frequency range 
about the pole. However, we will see in the next subsection that the sensitivity to g^T scale physics cancels 
when one sums over all ultra-soft gluon exchanges. For this subsection, where we are ignoring ultra-soft 
gauge bosons, one should also ignore g'^T scale contributions to the damping rate, and regard Fp as defined 
by an infrared cutoff at a scale /i such that g^T ^ <C gT. 
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K+P+Q 




P+Q 

FIG. 13. An individual loop in a ladder diagram with uncrossed rungs, with r, a vertex assignments 
consistent with a leading-order pinching-pole contribution. 




P 

FIG. 14. An individual loop in a cross-rung ladder. The light arrows inside the crossed loop illustrate 
the flow of frequency around the loop; the fact that frequency flows from a to r on both scalar propagators 
implies the absence of any "pinching-pole" enhancement. 



loop as shown. If both gauge boson propagators carry soft momentum, Q ~ ~ gT, then 
this is the scale of variation in frequency of the gauge boson propagators. And as before, 
this implies that over the frequency interval of width g'^T in which the P+Q and K+P+Q 
scalar propagators have nearly pinching poles, the gauge boson propagators are constant up 
to 0{g) corrections. Therefore they factor out of the frequency integral, which as we have 
already seen contains a pair of poles spaced 0{g^T) apart. So, just as for the end-most loop 
examined earher, we have an integral of form / dq^ GAdv{K+P+Q) GnetiP+Q), which is 
0{l/g'^T^). Therefore, every pair of scalar lines between uncrossed gauge boson rungs gives 
rise to a pinched frequency integral which contributes a l/g"^, as required to make uncrossed 
ladders contribute at leading order. 

Next, consider the case where two gluonic lines cross, as illustrated in Fig. Define uo 
to be the difference between the frequencies on the Qi and Q2 gluonic lines, uj = — q^i 
and (f as the corresponding average, (f = + \q2, so that dq^ dq2 = d(f duj. The a, r 
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FIG. 15. A soft vertex correction, which is also sub-leading. 



index assignments are already forced on us by the need to have pinching poles in the P and 
P+K scalar lines leading into this part of the diagram, plus rr assignments on the gluon 
rungs in order to obtain Bose enhancement factors on the gluon lines. The relevant part of 
the frequency integration thus becomes 

J doo GariK + P + Qi) Gra{P + Q2) = J doo G ^^,{K + P + Ql) Gnet{P + Q2) (4.11) 



/ 



duj 



We have written the denominator so that uo enters with positive sign in each term. Because 
uj is accompanied by —iV in both terms, the pole positions in the uj plane are all on the 
same side of the real axis, and the contour can be deformed to avoid all poles by a distance 
parametrically large compared to F. Hence there is no pinch contribution, and the crossed 
rung contribution is suppressed. An easy way to see that this will happen is to note how the 



frequency flows around the diagram, as illustrated in Fig. |T^. The frequency flows through 
both quark propagators from left to right, that is, from a to r; therefore (when one rewrites 
a retarded propagator with negative frequency as an advanced propagator with positive 
frequency) one gets a product of two advanced propagators, rather than one retarded and 
one advanced. 

Physically, the absence of crossed rungs can be understood most easily by Fourier trans- 
forming from frequency to time. The a, r assignments we have found necessary imply that 
the quark propagators, from one current insertion to the other, go between vertices in de- 



scending time order. Hence a diagram like Fig. 13 or Fig. TA, read from left to right, goes 
from later to earlier time. Two neighboring, uncrossed gauge propagators are time ordered; 
both vertices of the first propagator occur at a later time than either vertex of the second. 
On the other hand, when two gauge boson lines cross, the top vertices are in one time order- 
ing and the bottom vertices are in another. Since the time scale for a soft scattering event 
is ~ i/gT, while the typical time between scatterings is ~ 1/g'^T, it is natural that such 
overlapping scattering events are suppressed. 

With this experience it is now easy to see why diagrams containing a soft vertex correction 



like that in Fig. 15 are also suppressed. In terms of the added frequency integration variable. 
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both new quark propagators are traversed from a to r; the frequency integration will have 
two advanced propagators and there will not be a pinching pole. 

Hence, when all gauge boson lines carry soft 0{gT) momentum, ladder diagrams of the 
form shown in Fig. H contribute to the leading order emission rate, but crossed ladders of 
other forms or vertex corrections do not. 

E. Momentum exchanges softer than gT 

Our next task is to analyze the effects of gauge boson interactions carrying momenta 
which are parametrically small compared to the "soft" gT scale. We will proceed in two 
stages, first examining "very-soft" momentum scales which are intermediate between gT 
and g'^T, that is |q| = 0{g^T) for 1 < z/ < 2. We will argue that the previous analysis 
showing the dominance of uncrossed ladder graphs continues to apply, and then show that 
the net effect of these exchanges, at leading order, is merely to cancel the dependence on 
the contribution of momenta softer than gT to the quark self-energy. Then we will extend 
the discussion to "ultra-soft" g^T scale interactions, and argue that a similar cancellation 
continues to apply. 

To begin, it will be useful to review properties of the gauge boson propagator for momenta 
and frequencies which are in the range g^T <C |q| -C gT. In this domain, the longitudinal 
part of the gauge field propagator is screened at all frequencies, so that for any spacelike Q, 

G'^.,,^{Q)-llml = 0{llg'T'). (4.12) 

However, the transverse self-energy |31|] has frequency dependence which will be essential. 
For spacelike momenta with |q| <C (7^,0 

^^^Z.,^{Q) - {gT? , imHret,T(Q) - ^ {gT? . (4.13) 

If g° ~ g then the transverse self-energy is 0{g'^T'^) and the transverse retarded propagator 
is 0{l/g'^T'^), just like the longitudinal part. But if g° ~ |q|^/(5'T)^, then the self-energy 
is pure imaginary [up to corrections suppressed by q^/{gTy] and of order q^, leading to a 
retarded propagator which is of order 1/q^ — larger than the longitudinal propagator by 
a factor of {gTy/q^. This low (or "quasi-static") frequency regime dominates the spectral 
weight. In the very soft region, the equal time, transverse gauge field correlation function is 
given by its free theory value, 

z I £x e-'^- (A^(x, t=0)A'^(0, t=0)) = I ^ G^:{q\ q) ^ i (q) ^ , (4-14) 



^^For timelike momenta, there is a value of g° where the (real part of the) inverse propagator has a zero, 
leading to a peak in G^^^{Q)] this is the plasmon. However, the plasmon pole occurs at ~ mo = 0{gT) 
even for |q| <c; gT. For |q| ^ gT, the spectral weight carried by the plasmon pole is small. Moreover, 
the frequency range where the scalar propagators are all near their pinching poles always involves spacelike 
four- momenta on all giuon lines, so the plasmon pole plays no role at leading order in g for either soft or 
very-soft gluons. 
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with Pr^^(q) denoting a transverse projector, up to relative 0{c^ / g'^T'^) corrections from 
the longitudinal contribution and 0{g'^T /\c^) corrections from loops.|^ The loop corrections 
become 0(1) at the scale (7^T, reflecting the fact that physics at this scale is non-perturbative. 

When |q| = 0{g'^T), with u > 1, the d^q phase space will contribute a factor of (g^T)^, 
which is now parametrically smaller than the previous {gT^ soft phase space. We will see, 
momentarily, that only the low frequency, transverse part of the gauge boson propagator 
is large enough to produce a contribution in this momentum region which will affect the 
leading order emission rate. The longitudinal part of the propagator may be completely 
neglected. 

In the previous section, we integrated over loop momenta by doing the frequency in- 
tegration first, and argued that the only important frequency dependence [in the 0{g^T) 
interval which dominated the integral] was due to the scalar propagators. But as just noted, 
for very-soft momenta the gauge field propagator has significant frequency dependence on a 
scale of |q|^/ (gT)'^, which is of the order (or smaller) than g'^T for u > 4/3. Hence, we must 
modify the previous pinching-pole analysis to handle this regime where the gauge boson 
propagator depends strongly on To do so, first recall that the poles which nearly pinch 
in the scalar propagators GAdv{,K+P+Q) and G-Ret{P+Q) [cf. Fig. lie at = 
up to 0{g'^T) corrections. Instead of using gy and g° as independent variables in the loop 
integration, we will use q~ = q'^—q\\ and g". The essential point is that the scalar propagators 
depend only weakly on g° if q~ is held fixed.0 Consequently, the integral over at fixed 
q~, will be controlled by the frequency dependence of the gauge field propagator Glf!/{Q), 
not the scalar propagators, and receive its dominant contribution from the low frequency 
regime where q^ = 0{\q\^/g'^T'^). The nearly pinching poles in the scalar propagators will 
then appear in the q~ integral. 

The resulting parametric estimate for the relative size of a very-soft (|q| ~ g^T) gauge 
boson exchange is 



^'^The abse nce o f a mass term in the transverse part is, of course, a consequence of gauge invariance. The 
form of Eq. ( 4.14|) also has an interpretation in terms of dimensional reduction : the static limit of the 
transverse gauge boson correlator is given by "dimensionally reduced" 3-D Yang^^TvIills theory. 

^^For fixed q~ , a variation in by Sq^ changes the retarded propagator GRet{P+Q) by an amount 
SGnct{P+Q) - 25q° ^ GnaiP+Q) 

= 2<5qO [Gn,t{P+Q)? [ip-+q-){p++q+) - (p±+q±)' " ^Ret{P+Q)] 

= 2dq'> [Gnet{P+Q)? [{p-+q-) - {d/dq+) SRet(P+Q)] , (4.15) 

where q'^ = q^ + gy, and d/dq'^ = ■^{d/dq^ + d/dq\\). The momentum derivative of the self-energy is order 
g^T provided p = 0{T). And {p^ +q~) G-Rjet^P+Q) is order 1/T both within the pinching-pole region 
[where {p^+q^) is 0{g^T) and GRot(P+Q) is 0{l/g^T'^% and outside this region. Hence, for fixed q , the 
relative variation in the retarded propagator 5G-R,ct / G-Rct — 0{5q^ /T). Exactly the same estimate holds for 
the advanced propagator GAdv{K+P+Q). Provided q° varies by 0{gT) or less, this implies that if we treat 
the scalar propagators as depending on q^ and q\\ only through the combination g^, then we will be making 
at most an 0{g) error which is irrelevant to our leading-order analysis. 
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K+P 



K+P 





P+Q 




P+Q 



FIG. 16. Three contributions involving the addition of a very-soft gauge boson, |q| ^ gT, which cancel 
to leading order. 



q^T^ f 



GAdv(^+P+g)G'Ret(P+Q) 



27r 



2^ 



^^T^ X X (1//T^) X {T/q' 
1. 



(4.16) 



The g° integral, more explicitly, is /((ig°/27r)G^^(g°, q'||=q'°— g", q±) and coincides, to leading 
order, with the static result ( [4.14| ) for |q| <^ gT since the dominant contribution arises from 
0{\q\^/g'^T'^) frequencies which are small compared to |q|. Hence the difference between 
holding g|| or q~ fixed when integrating over g° is subdominant. The q~ integral over the 
product of scalar propagators gives an 0{l/g'^T^) result whose explicit form is just the 
appropriate adaptation of Eq. (|4.8|), namely 



J Zn 



4/ (p^+k) (r + i 6E) 



(4.17) 



[with 6E given in Eq. ( [4.9|) ], provided |q| <^ |p| and ~ p\\. 

The estimate ( 4.16| ) implies that a very-soft gauge boson exchange with momentum 
g'^T ^ |q| ^ gT can be just as important as a soft 0{gT) exchange. Of course, if one 
very-soft gauge boson exchange can make an 0(1) contribution to the emission rate, so 
will uncrossed ladders containing any number of such exchanges. Crossed rungs remain 
negligible in this very-soft regime for the same reason as before; if both propagators in the 
q~ integration are retarded or both are advanced, there is no pinching contribution. 

It is important to note that the very soft momentum region also makes an 0(1) contri- 
bution to the value of the on-shell scalar self-energy, which controls the size of the pinching 
pole contribution ([4.17|) . It therefore appears that momenta softer than gT are important 
at two points in the calculation: self-energies and uncrossed ladders. However, we now show 
that these two contributions from momenta q -C gT actually cancel to leading order in 
q'/g'T^ 

To understand this cancellation, consider the sum of the three sub-diagrams shown in 



Fig. 16, representing the contribution of a particular momentum q to the scalar self-energies 
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of a pinching pair of propagators and to a possible cross-rung between them. At this point, 
one should imagine having introduced a separation scale fi intermediate between g'^T and 
gT, and using propagators in which self-energy contributions from all scales larger than n 
have been resummed, but where self-energy corrections from scales softer than are treated 
as explicit loop corrections. 

To start, note that the momentum flowing along either leg on the right side of the 
fragments shown in Fig. |T^ must differ by iiQ for the case of very-soft gluon exchange as 
compared to either of the self-energy insertions. But this is a small correction to the rest of 
the diagram; adding a very-soft momentum q to p changes the parallel (to k) component of p 
negligibly since p\\ = 0{T), and changes the orthogonal component p^ ~ gT by a fractional 
amount which is at most 0{\q\/gT). In fact, once one averages over the direction of q_L, 
the relative size of surviving corrections is only 0(q^/p^) ~ / g'^T'^. Hence we need only 
check whether the contributions of the sub-diagrams pictured cancel. These contributions, 
including the and g° frequency integrals but neglecting common overall factors, are 

X \^{2K+2P+Q)^{2K+2P+Q),[G^A,{K+P)fGKAAK+P+Q)GKct{P) 
+ {2P+Q),{2P+Q),Gp^UK+P) [GRet(P)]'GRet(P+Q) 

+ {2K+2P+Q)^{2P+Q),G^UK+P)G^UK+P+Q)G^,,,{P)G^,,{P+Q)^ . (4.18) 

At leading order, all the terms of form {2K+2P)^ are parallel with K and may be replaced 
with their magnitudes times K'^, where K = (l,k). Further, |q + k + p| |k + p| ^ + 
The scalar propagators depend weakly on g° for fixed q~, so the g° integral may be performed 
as discussed above. Next, we may approximate each scalar propagator by just the relevant 
nearly-pinching single pole. For convenience, we will also define a shifted frequency variable, 
p/o ^ pO _ sign(p||), and then drop the prime so that 

^"""'^^^ = - [po + ,r/2]2 ^ ~2EpSign(p||) po + ,r/2 ' ^^'^^^ 

1 11 
Ga^AK+P) = ■^2^^Z:j^^^—^ ^ ~2Ep+ksign(p||+A;) + SE - lY 12 ' ^^'^^^ 

and likewise for all the other propagators. The key observation is that the q^ depen- 
dence of the difference between the pole positions of G'Rct(-P) and G-^^^iP^Q), or between 
G'Adv(-^+-P) and (jAdv(-^+-P+<5)5 are negligible provided q <^ gT . In other words, 

bE = sign(]9||) - sign(p||+/i;)Ep+k + 

- sign(pii) Ep - sign(p||+g||+/i;) ^p+k+q + + q\\ , 

- sign(p||+g||) Ep+q - sign(p||+A;) Ep+k + A;° - q\\ , 

- sign(p||+g||) ^p+q - sign(p||+g||+fc) ^p+k+q + . (4.21) 

This amounts to neglecting the contribution of q^^ + 2pj_ • qj_ to (p±+q±)^. The corrections, 
relative to the size of bE^ are of order q^/(7^T^ after averaging over the direction of q_L, 
since p^^ ~ 9^T'^ and p_L ■ q_L vanishes upon angular integration. 
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Under these approximations, and recalling that <^ q ^ q\\ for very-soft q, expression 
( 4.181 ) becomes 



1 



{p^+iT /2){p^-iT /2+6E) 



+ 



{p^-iV/2+5E){p^-q\\-iV/2+5E) {p^+iV /2) {p^-q\\+iT /2) 
1 1 



{p^-q\\-iTl2+6E) {p^-q\\+iT/2) 



(4.22) 



We now perform the p° integral. The first two terms have one pole on one side of the contour 
and three on the other; the last term has two poles on each side, so it gives one contribution 
for each residue. The result is proportional to 



1 



1 



iT-5E (iT-6E 



1 

+ — 



{ir-6E) q 



+ 



1 



'iV-SE 



1 



{{r-6E) 



1 



0. (4.23) 



Corrections to this cancellation arise from the corrections to the approximate equalities 
( [4.21| ), as well as corrections from neglecting Q (relative to P) in other parts of the diagram. 
All such corrections are suppressed, after averaging over the direction of q^, by at least 
0[q^/(5'T)^]. In particular this cancellation is not sensitive to corrections to the static cor- 
relator ([4.14|) , and so is unaffected by the loss in perturbative calculability of this correlator 
when |q| = 0{g^T). 



To be complete, we should also check that the three diagrams of Fig. ^ have the same 
group theoretic factor in a non-Abelian theory. This turns out to be trivial. In any diagram 
where the gluonic lines never cross (in the sense we have used above), so called "rainbow 
topologies," the non-Abelian group factor isQ (Cp)*^'"""*^. 

The physical interpretation of the cancellation ( [4.23[ ) is as follows. There is an 0(1) 
chance that during the order 1/g'^T time required for the photon emission process to com- 
plete, there will be a scattering with exchange of g -C gT momentum. However, such a 
scattering disturbs the quark (or scalar) too little to disrupt the emission process. The 
self-energy correction accounts for the removal of the quark (or scalar) from the initial mo- 
mentum state, and the ladder rung accounts for its contribution in the new scattered state. 
This interpretation will become more manifest when we resum ladder diagrams in the next 
section. In a relativistic context, this cancellation is the same as that discussed by Lebe- 
dev and Smilga |^ in the context of electrical conductivity (see also Ref. |2^), but the 
basic observation goes back decades.0 However, these works address a different kinematic 



^^Note that this would be false, and the cancellation ( 4.23| ) would not occur, if we were considering color 
conductivity 

^''A single small-angle scattering does not substantially affect transport processes such as electrical conduc- 
tivity. In the non-relativistic context, this insensitivity to small angle scattering appears as a suppression 
factor of (1 — cos 6*) in what is sometimes called the "transport cross-section" p^]. This factor is relatively 
easily derived from the Boltzmann equation. Diagrammatically, it emerges as a cancellation, when sumrning 
ladder diagrams, between rung and self-energy insertions. An early discussion may be found in Ref. [pSj ; 
see also Sec. 39.2 of Ref. M. 
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Q 




FIG. 17. Two ways of adding a g^T line, denoted by the thin gluon line with momentum Q, which will 
cancel at leading order. 



regime than the one relevant here (they consider extremely soft /c, whereas we have hard 
but lightlike fc), and the range of exchange momentum for which the cancellation occurs in 
the current context is not the same as in these references. In particular, in our application 
an exchange momentum |q| ~ gT is enough to disturb the collinearity condition which must 
be satisfied for the poles to pinch. This is not the case for softer k. 

Having found that the leading order photon emission rate is insensitive to the effects 
of very-soft {g'^T <^ |q| <^ gT) gauge bosons, one should naturally expect the emission 
rate to be equally insensitive to the ultra-soft g^T scale. (Otherwise, one would have a 
peculiar situation in which a physical quantity, at leading order, depends on the dynamics 
of two different momentum scales, gT and g'^T, but not on intervening scales.) Certainly, 
the physical picture just sketched does not suddenly change when q reaches the g^T scale. 

However, some readers may want a more deductive demonstration that this insensitivity 
remains valid at the g^T scale. To this end, first observe that nothing changes in the 
analysis of the diagrams of Fig. |16| if the momentum q is in the ultrasoft 0{g'^T) regime. 
The cancellation among these diagrams remains applicable at the non-perturbative g^T 
scale. However, in this regime, the previous arguments which ruled out crossed diagrams 
and vertex corrections no longer apply; one should verify that a similar cancellation also 
occurs between diagrams of the form shown in Fig. The demonstration is very similar 
to that for uncrossed graphs. 

It is sufficient to show that when we sum over the two possibilities of attaching the 
left vertex of the g^T gluon line to either the upper or lower quark lines, there is an 0{g) 
cancellation. The same argument will apply, independently, for the attachment point of the 
right vertex, leading to an overall cancellation by 0{g'^). The cancellation we will find also 
applies if the other end of the g^T line attaches to another ultrasoft gluon line, so it will 
rule out more complicated ultrasoft corrections as well. 

First note that the addition of the ultrasoft q ~ g^T line does not modify, at leading 
order, any of the integrations over loops which the ultrasoft propagator skips over. This 
is because the transverse momentum it carries is smaller by g than the typical transverse 
momentum already in the loop, and the longitudinal momentum is smaller by g^. Therefore 
we need only focus on the pair of pinching poles which are augmented by one more propagator 
when (one end of) the ultrasoft line attaches. The sum of the two contributions at fixed q, 
after performing the integration and suppressing factors of g and group factors, is 
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FIG. 18. A diagram fragment in which an ultra-soft gluon line connects two soft gluon rungs. Such 
diagrams are suppressed by 0{g^). 



+ 2P^ Ga^AP+K) Gnet{P-Q) G'Ret(P) 



4:p\\{k+} 



q 



+ 



1 



T 
q2 



{p^-iT/2+SE) {p^-iT/2+6E+q\^) {p^+iT/2) 
1 

{p^-iT/2+6E) (p0+ir/2-g||) {p^+tT/2) 
1 1 



(tT-dE) {tr-6E-qii) {{r-6E) {tT-dE-qii 



(4.24) 



and again the two terms cancel. (In the first step, we shifted p° by Ep and retained only 
the pinching pole terms, exactly as before.) And just like before, one can verify that non- 
Abelian group factors do not affect this cancellation (because the currents at either end of 
the full diagram are gauge invariant). 

It is elementary to see that there are also no leading order contributions when an ultrasoft 
line attaches to one of the soft exchange gluons. This is because the vertex involves a 
derivative coupling and only an order gT momentum runs through the soft gauge boson 
line. To be definite, consider an ultrasoft gauge boson connecting two soft gauge boson 
lines, as shown in Fig. ^3. No new pinching poles are added. There is an explicit from 
the new vertices, a g^ from the derivatives at the vertices, a g^ from the ultrasoft phase 
space, a.l/g^ from the frequency integration over the ultrasoft rr propagator, and l/g"^ from 
two new ra soft gauge boson propagators. Putting it together, there is a g^ suppression. 
Physically, this reflects the fact that the soft exchange gluons only propagate a distance 
1/gT, much shorter than their mean free path to interact with the ultrasoft field which is 
of order 1/g'^T. 

Therefore, although scales softer than gT appear to contribute if one examines individual 
diagrams, when appropriate combinations of diagrams are summed, contributions from these 
scales cancel. There is therefore no danger in temporarily imposing an infrared regulator, 
whose form we need not specify, which eliminates the sensitivity of the decay width Pp 
to ultrasoft physics and allows the manipulations of the next section. After summing the 
required diagrams, the infrared regulator may be removed without re-introducing any poor 
behavior [i.e., infrared sensitivity). 
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F(P,K) 



M(P,K,Q) 



FIG. 20. Correspondence between graphical elements and the symbols appearing in the integral equation 
which represents the sum of all ladder graphs. 



V. SUMMING LADDER DIAGRAMS 

Summing the contributions of all ladder diagrams requires a procedure analogous to 
the formulation of a Bethe-Salpeter equation; one cuts off the photon vertex from one end 
of every diagram, recognizes that the resulting sum is just a geometric series, and writes 
down the linear integral equation whose iteration generates this series. This is illustrated in 
Figs, m and EUl 



To carry this out explicitly, we first introduce some notation, chosen to mimic the treat- 



ment of Jeon 1 40 1 as much as possible. We define a bare photon vertex X^{P] K), as well as 



its transverse projection X^(P; K), 

I^'iP; K) = {2P+KY , X^(P; K) = {2P+KY^ = 2P{^ , (5.1) 

which represent the vertex factors associated with the attachment of an external photon. The 
full X'^ will be used to find the complete current- current correlator, while Xj^ will generate 
the transverse projected correlator, which turns out to provide the most economical form 
for extracting the photon production rate. We will want both in what follows, in order to 
check explicitly that our leading-order result for the current-current correlator is transverse. 

We introduce J-'{P; K) to denote the pinching-pole part of a pair of scalar propagators, 
with r, a indices appropriate for evaluating Gaarr] 



pinching— pole 
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-1 



4p|| {p\\+k) T + i6E 



An 6 



2/ + A;°-Epsign(p||) -Ep+kSign(p||+A;) , (5.2) 



where the second expression is the pinching pole approximation, and F is shorthand for 
|(rp+rp+k). The delta function represents the product of nearly pinching single poles, as 



in Eq. (O), and so should really be smeared out over the 0{g'^T) width of the pinching-pole 



region. But since the exchange momenta are 0{gT), approximating the product of (off-axis) 
poles by the indicated delta function is sufficient. We have also, once again, used the fact 
that p\\ will be 0(T) while is 0{gT) to write p\\{p\\+k) in the prefactor instead of EpEp^^ 
(times an appropriate ± sign). 

Next, we introduce a linear operator A4 which will represent the addition of a gauge 
boson cross-rung with rr indices. The kernel of this operator is 

MiP, Q; K) = iCf 2P^ 2{K+P), ^^^^(Q) . (5.3) 



The overall factor of i reflects the i's associated with each vertex (|3.11| ) together with one 



—i for the propagator to compensate for the i included in the definition ( |3.4|) of We 
have already made a. Q <^ P, K approximation in writing 2P^ rather than {2P+Q)^, etc. 
We may also use the lightlike nature of K and (near) collinearity of P and K to further 
simplify the kernel to the form 

M{P, Q; K) = zg'Cp ipnipn+k) K^K, G^^::{Q) , (5.4) 

with = (1, k) as before. Note that the group factors are trivial since none of the gluon 
lines cross; there is just a factor of g^Cp for each gauge boson line. 

Finally, we introduce an effective vertex T>^[P] K), which represents the sum of all num- 
bers of cross-rungs and pinching-pole propagators placed in front of the photon vertex X^, 

= P' + MJ^I^" + MJ^MJ^V + MTMTMTX^ + • • • . (5.5) 

More explicitly, 

V^{P;K) =I^{P] K) 

+ / M{P,Qi;K) J^{P+Q^-K)I''{P+Q^-K) 

JQi 

+ f M{P,Qi;K)T{P+Qi;K) [ MiP+Qi,Q2;K)J^{P+Qi+Q2;K) 
Jqi Jq2 

X j^(p+gi+g2; K) 

+ ■■•. (5.6) 



As illustrated in Fig. this is alternately expressed as the linear integral equation 



P^(P; K) = P'(P; K)+ f MiP, Q; K) T{P+Q- K) V^'iP+Q; K) . (5.7) 

JQ 

As also illustrated in Fig. the resulting LPM contributions to the full current-current 
correlator are proportional to the real part of an "inner product" of XJF with P, 
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WHUK) = -J ^ 2n{p'+k') [l+n{p')] V\P- K) Re [^(P; K)V\P- K)] . (5.8) 

The overall minus sign is present because the population functions, as given by Eq. ( [4.5| ), are 
n{jp +kp)n{—p^) , and n{—p^) = — [l+n{p^)]. The fact that only the real part enters follows 
from the discussion after Eq. ( |4.5| ). Although it is not obvious from its structure, Eq. ( ^.8|) 
does yield a current- current correlator W^^lK) which is symmetric under interchange of n 
and u. This will be shown in Sec. IVl B. 



As it stands, Eq. (|5.71 ) is a four- dimensional integral equation characterizing the P depen- 
dence of V^'iP; K). However, the delta-function in J'(P; K) implies that J'(P; K) Vf{P; K) 
has support only for a single frequency for a given spatial momenta p (and K). Conse- 
quently, it is convenient to define 

r (p; k) ^ -4p|| (pii+A;) j ^ ^(P; K) V^{P; K) , (5.9) 
in terms of which the current- current correlator (15. 81) becomes 



W^UK) = (2p.+/,.)Rer(p;k). (5.10) 

J (277)^ 4p|| (fc+pii) 

In this equation, and in the equations which follow, the time components of 4- vectors are to 
be understood as dependent parameters, fixed by the locations of pinching pole pairs; so in 
the above, 2p^+kP = ii^p sign(p||) + £'p+k sign(p|j+/c). 

The integral equation ( |5.7|) may be reduced to an equation just involving f^{p; k). To 
do so, multiply Eq. ( |5.7D by J^{P]K), integrate over and then multiply both sides by 
4:pi\{p\\+k){T+i6E). This yields 

2P^+K^ = {z6E + r) rip; k) - I ^ C(q; k) ^ (p+q^; k) , (5.11) 

where 

C(q;k) =/Cf /^27r<5(g°-g||) [-z G'^f;(g)i^^i^J . (5.12) 

The replacement of (p+q) with (p+q±) in the first argument of inside the ci^q integral 
is justified at leading order because q\\ is negligible compared to p\\, and the variation of 
/'^(p; k) with p|j occurs on a scale of T. [In contrast, the /'^(p;k) varies with p_L on the 
gT scale.] We have also made the approximation E^^^—E^ = gy sign(p||), valid at leading 
order, to turn the delta function in J^{P+Q, K) into S{q^—q\\). 

Next, we need to insert the actual form of the damping rate F. It is related to the 
imaginary part of the on-shell retarded self-energy, 

p° Tp = /Cf Im J 4P^P, [-zG^^;(Q)] G^,,{P+Q) . (5.13) 

Everything is real except the retarded propagator; its imaginary part is half the quark 
spectral weight. 
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27r(5(g°-g||: 
4pO 



(5.14) 



where we used the on-shell condition p° = _E'pSign(p||) and di^p+q/rfq ~ p ~ k, dropped the 
irrelevant pole at ^ — -E'pSign(pj|), and took the quasiparticle width to be small compared 
to q, which is legitimate for q ~ gT. (Remember that the g'^T scale is irrelevant and we 
have implicitly regulated it.) Consequently, 



(5.15) 



and therefore we can rewrite Eq. ( |5.11D as 



2P^+ir^ = «5E/^(p;k) + 



(2vr) = 



C(q;k) r(p;k)-r(p+q^;k) 



(5.16) 



This is structurally similar to a Boltzmann equation; in real space the i6E term can be 
viewed as a gradient (flow) term, and C(q; k) as the rate for scattering with momentum 
transfer q. The self-energy has turned into a loss term, accounting for the scattering of 
excitations out of momentum p, and the ladder exchanges have generated a compensating 
gain term, accounting for scattering into momentum p from momentum p+q. 



It remains to present the explicit form of the "rung" —i Gt^!;!{Q)K^K,y 



This contraction 



is gauge invariant, when multiplied by S{q'^—qi\) [as in Eq. (|5.12|) ], and we are free to express 



in the most convenient gauge, which will be Coulomb gauge. With this choice. 



Ret 



iQ) 



- nL(Q) ' 



Ret 



iQ)=0: 



q2_(^0)2 + nT(Q) 



(5.17) 



For spacelike Q <^T [enforced by the delta function in Eq. ( |5.12|) ], 

nL(Q) 



nT(g) 



UIt 



m 



2q q — q' 







.vrg 



2g2 



q-q 



.Tcq 
~^''Tq 



(5.18) 
(5.19) 



with q= |q|, and mo the leading-order Debye mass. The overall signs in Eq. ( p.l7| ) reflect 
our metric convention; in particular, g^^ = — 1 appears in the numerator of G^^ . Recalling 
that is a null vector with unit time component, using Eq. ( p.4| ) and Eq. ( |3.10| ), and 
approximating nb{q^) ':^T /q'^ (valid at leading order for g° ~ gT)., one flnds 



iG':';{Q)k^k, 



2g 



[i-iqVqm-iq\\/qf 

(gO)2 + nT|2 



(5.20) 



This completes the derivation of Eq. (|1.6|) 
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VI. DISCUSSION 



A. Symmetry 



In thermal equilibrium, the Wightman current- current correlator ( ^.21 ) is symmetric un- 
der interchange of indices, W'^'^{K) = W^'^{K). This relation should follow automatically 
given the structure of the diagrams which we summed, but it is not manifest in our result 
( ^■10| ) for the correlator. Hence, it is worthwhile to check that this symmetry nevertheless 
still holds. 

The integral equation ( 5. IE ) controls the dependence of f^{p] k) on p^, for fixed values 
of p|| and k. To write this equation more abstractly, consider the vector space of smooth 
functions of p_L equipped with the natural inner product 



{9\f) 



9{p±rf{p±). 



(6.1) 



One may define linear operators 6E and C in the obvious way so that the integral equation 
( p.l6| ) becomes 

\I'^)=\t6E + c]\n, (6.2) 
and Eq. ( p.lO| ) turns into an integral over with an integrand proportional to Re (X^|/'^), 



4p|| {k+i 



(6.3) 



[Explicitly (J'^) = 2P/^+ir^ C\f) = /rf3q/(2vr)3C(q;k)[/(px)-/(p^+q^)], and 6E is 
just a multiplication operator by 6E, all for fixed given values of p\\ and k.] 

The operators 5E and C are both real and symmetric, and hence Hermitian; for the 
collision operator C this is a consequence of the ~c[± symmetry (or rotational invari- 

ance) of the kernel ( |5.2CI| ). Moreover, (SE)"^ is positive definite on the space of normalizable 



functions, which implies that the operator 
write 



i6E + C 



is invertible. Consequently, one may 



ff") = i5E + C 11") . 



-1, 



If A4 denotes the real part of the inverse. 



M 



i6E + C 



+ 



i6E + C 



-1 



then 



and 



dp\\ 2ra(/+fc°) [l+ra(/)] 



2tc 4p|j (fc+P|| 



{i^\M\r). 



(6.4) 

(6.5) 
(6.6) 

(6.7) 



Because the operator Ai is, by construction, real and symmetric, this form shows that our 
integral equation does yield a current-current correlation function which is symmetric under 
interchange of indices. 
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B. Trans versality 



The exact current-current correlator is transverse, Kf^W^'^{K) = 0, and this imphes that 
the photon emission rate may be computed either by summing over all external polarization 
states with g^^W^i,{K), or by summing only over the two transverse physical polarization 
states with Y.a=v2^f\a)^^lv{K). 

Our result (|5.10|) is not manifestly transverse, and it is certainly appropriate to ask 
whether this property nevertheless holds. If we had exactly evaluated all of the diagrams 
we summed (and included all "rainbow" diagrams in the scalar self-energy), then our result 
would be exactly transverse. However, we used various approximations in evaluating these 
diagrams which, in the kinematic region of interest, were valid to leading-order but neglected 
effects suppressed by additional powers of g. [In particular, this occurred in the steps leading 
to Eq. ( |5.11|) and Eq. ( |5.12|) where q\\ dependence in /'^(p-|-q;k) and O(q^) shifts in the 
pole position of ^{P+Q] K) were neglected.] Consequently, it is quite possible that our 
leading-order approximation ( p.lO| ) to the correlator will not be exactly transverse. What is 



essential, however, is that any non-transverse part of W^^ be parametrically small relative 
to the transverse part. 



To examine this issue, contract both sides of the integral equation (|5.16|) with to 
obtain 

{K' + 2P,K^) = {tSE) K,r{p; k) + J ^ C(q) [K.rip; k) - TT^/^ (p+q^; k)] . (6.8) 
If one could neglect the collision term altogether, then the solution would trivially be 

ir^r(p;k) = ^Tr^ = ^ £;pSign(p||) + Ep+ksign(p|,+A;) , (6.9) 

where the second equality depends on the pinching-pole condition + = i?psign(p||) + 



-Ep+k sign(p||+/c) plus the exphcit form = yp^+m^ for the quasiparticle dispersion 
relation which neglects 0{g^T) corrections due to momentum dependence in the self-energy. 

Now, 

Epsign(pii) + ^p+kSign(p||+A;) = 2p\\ +k + 0{pl/p\\) , (6.10) 

and the leading 2p\\+k piece cancels exactly when inserted into the collision integral in 
Eq. ( |6.8| ). Consequently, the collision integral is a small perturbation relative to the other 
two terms in Eq. (|6.8| ). Iteratively correcting the trial solution (|6.9| ), one easily sees that 
X^/'^(p;k) differs from the trial solution only by terms of order g'^T. Because the trial 
solution (|6.9|) is pure imaginary, this means that the real part of K^/^(p;k) is suppressed 
relative to the imaginary part by 0{g^). And since the current-current correlator ( |5.10| ) 
only depends on the real part of /^(p; k), this means that the non-transverse part of (our 
approximation to) W^^{K) is at most 0{g^T^). Note that it was essential to this argument 
that the collision term has both a gain and a loss term, of equal size. 

In contrast, the transverse part of /'^(p;k) is of order p^/6E ~ 1/g [because both 6E 
and the collision operator C are 0{g'^T)], which leads to an 0(1) transverse contribution 
to W^'^{K). [The extra power of g comes from the factor of 2Pj^ = 0{gT) which appears 
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in the integrand of Eq. (|5.1CI| ) when one speciahzes to the transverse part.] Consequently, 
the longitudinal part of our result for the current- current correlator is suppressed by 0{g^) 
relative to the transverse part.0 



C. UV and IR behavior 



We next briefly examine the behavior of the integral equation ( |5.16|) in the limits of large 



and small transverse momentum p^, and large and small momentum transfer q, in order to 
confirm that the photon emission rate is finite and depends (at leading order) only on gT 
scale physics. Since we have just verified transversality of the correlator, it is sufficient to 
consider only the transverse components of /'^(p;k). For p_|_ = 0{gT) and p\\,k = 0{T), 
recall that both the collision operator C and 6E are 0{g^T). Therefore both the real and 
imaginary parts of are 0{l/g). Now if p_L ^ gT, then the iSE term in Eq. ( |5.16| ) 
dominates over C, leading to 



mpr, k) --2^-^ -4. . (6.11) 



Hence Im/^ vanishes as 1/|p_l| for large p_L. To determine the asymptotic behavior of the 
real part, we must expand in C/6E, leading to 

Re/^(p^;k)~-^|^C(q)Im[/^(p)-/^(p+q^)] . (6.12) 

At large q, the collision kernel C(q) oc J^o G!^!;! {Q)\g^^=gO ~ This renders negligible the 

contribution from |q^| ^ |p±|- The contribution to Re from |q_L| ~ |p^| is of order 
p±/{p^±)- To estimate the contribution from |q_L| <^ |p_l| we may make the approximation 
that, on angular averaging, 

fl(p) - /r(p+qJl..,,„„,,„„,, ~ -il fiip) . (6.13) 

to conclude that the contribution from |q^| <^ |p_|_| is at worst logarithmically sensitive to 
small q (and cut off at |q^| ~ gT), so the small qj_ contribution changes the |qj_| ~ |pj_| 
estimate by at most a logarithm. Hence, ignoring logs, 

Re/^(px;k)~^ ^^-^ . (6.14) 



This decreases sufficiently rapidly with p_L to ensure that the d'^p± integration in Eq. (p.lO|) 



converges at large p_L and is dominantly sensitive to p_L = 0{gT). One may easily see that 



^®In fact, our result (5.10) for the current-current correlator is exactly transverse when the photon mo- 
mentum is precisely on-shell, = 0, if one uses a quasiparticle dispersion relation containing just 



the first two terms in the large momentum expansion, E-p 



\P 



i(pi+m2 



)/\P 



In this case, 



i4r^/'^(p; k) — 2P^K^Ji5E is identically equal to i{2p\\+k), with no sub-leading corrections. Hence the 
collision term in Eq. ( p^ ) vanishes identically, as does the real part of K^f^^p; k). 
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there is no problem at small p_L, because both 5E and the collision term are regular at 
= 0. In fact, /j^ vanishes linearly in at small p^. 

The large q behavior of the gauge field correlator, GrriQ) ~ q^^, ensures that the 
collision term is well behaved at large momentum transfer. But the small momentum transfer 
behavior requires a bit more care. First note that, for all p_L, the value of /j^(p;k) varies 
with p± on a scale of order gT or larger, and so is slowly varying on smaller scales. This 
follows because the collision term smooths on this scale, and 1/SE is everywhere smooth 
on this scale. [Note that 1/6E is smooth even when p±_ <^ gT due to the thermal mass 
~ g^T"^ contribution to SE, as seen in Eq. Therefore, for small momentum 



transfer q <^ gT, we may again approximate the difference /j^(p;k) — /j^(p+q_|_; k), after 
angular averaging, by a Laplacian, 

/^(p;k)-/^(p+q^;k)] ~ Vj^ /^(p; k) . (6.15) 

J angle avg. 

Hence, this difference vanishes like g^. q_L — ^ 0. This is what will control the small 
q behavior. As discussed in subsection |1V K| , for q <ti gT, G^^(Q) is sharply peaked in 
frequency about g° ~ 0. At small q^, using Eq. ([4.14|) , one has 

^2.5(gO-gil) HG-(g)] ^ ^^1^ . (6.16) 
(27r)^ q_L 

Integrating this over q^ would generate a logarithmic IR divergence (or rather, logarithmic 
sensitivity to g^T scale physics), exactly as found in the damping rate F. However, our 
collision integrand also contains the difference ( |6.15D so the relevant small q^ behavior is 



i 



d^q^ 1 

q2«g2T2 (27r)2 



qIK. /±(P5 k)] ~ mp; k) / dq^ . (6.17) 



This is IR safe by two powers of q, and confirms the irrelevance of momentum transfers 
which are much smaller than 0{gT). As expected, it is a cancellation between the gain and 
loss terms, which arise from the gluonic cross-rung and self-energy, respectively, which leads 
to this good IR behavior. A treatment which includes only some of the relevant diagrams 
would fail to find this cancellation. For instance, including the self-energy but failing to 
sum ladders would miss the gain term, leaving logarithmic sensitivity to the g^T scale. 
This explains the conclusions of (Note that such a treatment also fails to produce a 

transverse current-current correlator.) 

Finally we remark on the small py behavior. In this regime, the scalar result ( p.3| ) 
behaves as J{dp\\/p\\^) Re (p_L | /_L(p;k)), with one factor of l/p\\ explicit in Eq. (|6.3|) and 
one coming from the Bose distribution nb{p^) ~ T/p\\ for near on-shell P. This will lead 
to a small p\\ divergence (or more properly, linear sensitivity to the inverse thermal mass) 
unless Re /^(p; k) vanishes as p\\ — > 0. This is the case, because 6E ~ l/p|| [cf. Eq. ( |1.13| )], 
which becomes large for small py. When 6E becomes large, our previous analysis leading 
to Eq. ( |6.14| ) shows that Re /j^(p; k) oc py^, which renders the p\\ integral in Eq. ( p.3[ ) well 



behaved. The analysis for fc-t-py <C T is exactly the same. For the case of fermions instead 
of scalars, the expression in Eq. ( |1.3| ) has one further explicit inverse power of py, but the 
Fermi statistical function does not have singular small energy behavior, so the overall small 
py (or small p\\+k) behavior is the same. 
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VII. EXTENSIONS 



A. Fermions 

It remains to extend our treatment to the physically relevant case of quarks; that is, 
charged, nearly massless fermions instead of charged scalars. The entire power counting 
discussion is essentially unchanged. The only differences arise from the Dirac structure of 
the massless fermion propagator which has the form 

^ (7.1) 



with a thermal self-energy ^(P) which does not respect Lorentz invariance (because the 
plasma defines a preferred frame). At large momentum, |p| ^ gT, 



2 



[Re S,(P)]^ = -[Re E,{P)Y + [Re E(P)]^ = -f^ , (7.2) 

with = ^g'^CpT^, the same asymptotic mass as in the scalar case. Therefore the mass 
shell condition (P + Re E)^ = is satisfied for a timelike 4-momentum, P^ + = 0, and 
the dispersion relation is the same as for a scalar with vanishing vacuum mass. 

Because we are regarding the zero-temperature fermion mass as negligible, the theory is 
chirally invariant, and this chiral symmetry is unbroken in the high temperature plasma.^ 
This means that we may consider separately the contributions to photon emission from right 
and left handed two-component (Weyl) fermionic fields. For the right handed component, 
the 7 matrices are replaced by 2 x 2 matrices where ^ 1 and a* are the Pauli matrices; 
for the left handed components, the relevant matrices are which differ just by changing 
the sign of a^. It is sufficient to consider the right handed component, as the unpolarized 
emission rate from the left handed component is exactly the same. 

Rotational invariance guarantees that the spatial self-energy 5](P) is parallel to p. This 
allows us to decompose the fermionic propagator as 

1 _ i(aO + a-p) |(a°-a-p) 



a^^{P + S)^ |p+S| - (/+S0) |p+S| + (/+S0) ' 



(7.3) 

The numerators of Eq. (|7.3|) are helicity projection operators. It will be convenient to let 
u{p) and v{p) denote the positive and negative hehcity eigenspinors of p ■ a, respectively, 
normalized to 2|p|. Hence, 

IpI (a° + a-p) =M(p)M(p)t, IpI (a° - a ■ p) = t;(p)t;(p)t . (7.4) 



^^Alternatively, one may say that the self-energy does respect Lorentz invadapce but depends on the plasma 
4- velocity u'^ as well as the fermion momentum F^. For a discussion, see 

■^"Non-perturbative effects associated with the U{1)a axial anomaly are irrelevant for our perturbative 
analysis. 
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Only one of the two pieces of the propagator ( |7.3| ) will contribute to the pinching poles 
in our ladder diagrams, so at leading order the piece with a non-pinching pole can simply 
be dropped. To be specific, consider bremsstrahlung, so that p\\ > 0. For frequencies near 
the pinching pole position, each fermion propagator G'Adv(-P) may be approximated by 



2p|| (po - - zT/2) 



(7.5) 



for p|j > 0. This is exactly the same as the scalar propagator except for the factor u in 
the numerator. To obtain an expression as similar as possible to the scalar case, we may 
associate the u and spinor with the vertices on either side of the propagator, instead of 
the propagator itself. The gluon vertex is essentially unchanged, because 



(7.6) 



with p° = |p|, up to 0(|q|) corrections which may be dropped. For the external photon 
vertex we need only the transverse components. It is convenient to work in terms of the two 
circular polarizations, cr"*" = [a^+ia"^ )/V2 and a = {a^—ia^)/\^, for which 



u 



2(P||+fc)^ 
P\\ (PW+k) 



Pj 



u 



^(p+k) Or-M(p) 



P\\ (PW+k) 



pi 



(7.7) 



(Switching to the other helicity just interchanges the roles of and a~ .) The corresponding 
scalar quantity is just {2p±Y. 

In the fermionic case, each diagram has an overall statistical factor of n f{k+p)n f{—p) = 
nf{k+p) [l—nf{p)], replacing the factors of nh{k+p) [l+ni,{p)] in the scalar case. Hence, the 
leading order contribution of a right (or left) handed fermion to the photon emission rate, 
summed over polarizations, differs from the scalar contribution by a ratio of 



fermion 
scalar 



2p|| ip\\+k) 



r^fjk+Pw) [l-Uf 
Mk+P\\) [l+nb{p\\)] 



(7.8) 



A Dirac fermion has two chiralities so there is a further factor of 2. Exactly the same 
result is obtained in the case of inelastic pair annihilation. (Note that for inelastic pair 
annihilation, the sign of the first factor is negative, as —k < p < 0; but this compensates 
for a sign difference in the population functions, because [l+n{,(]9)] = —nb{—p), whereas 
[1— n/(p)] = +nf{—p).) This completes the derivation of the fermionic result in Eq. (pTsj ). 



B. OfF-shell photon production 

For some purposes, off-shell photon production can also be of interest. For instance, the 
behavior of Wnu{K) for timelike K is important for dilepton production. Our analysis up to 
this point never assumed that the photon momentum is exactly null; for ^ one must 
merely include the appropriate photon frequency shift k^—k when evaluating the energy 
shift ( p.. 71 ). This has one noteworthy effect, however. With a null photon momentum, 
the energy shift 6E is strictly negative for pair annihilation [p||(p||+fc) < 0] and strictly 
positive for bremsstrahlung [p\\{p^\+k) > 0]. For a non-null photon momentum, k^—k must 
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be added to 5E and this implies that, for some range in py, the energy shift 5E can vanish 
for a particular value of p^^. This is the point where one or the other of the processes 
(bremsstrahlung or pair annihilation) is kinematically allowed without any accompanying 
soft scattering off the plasma. 

When g^T <^ \k^—k\ ^ T, this kinematically allowed point occurs when gT <^ \p±\ <^ 
T. In this regime, our treatment is valid but needlessly complicated. At leading order, 
the existence of soft scattering processes is unimportant, and one may drop the collision 
term altogether and merely apply an appropriate ie prescription to deal with the pole in 
f'^{p±]P\\, k) where 6E vanishes. This results in a standard 2 1 particle pair annihilation 

> k) or DIS (A;° < k) treatment. 

If \k^—k\ ~ T, then the value where SE vanishes occurs, for typical p\\, at ~ T^. 
Near-collinear processes (with or without associated soft scatterings) cease to dominate the 
emission rate. In this case, all our p^^ -C py^ approximations break down, and our treatment 
not only can be superseded by a simpler one, it must be. 

Finally, if |A;°— A;| < 0{g'^T), then the detailed structure of the collision terms remains 
relevant, and one must solve our full integral equation ( |1.6| ) to obtain the correct leading 
order emission rate. 

C. Softer photons 

Up to this point in our analysis, we have focused on the case of a hard photon momentum, 
/c ~ T. We now explore the modifications required when this condition is relaxed and the 
photon momentum becomes small compared to the temperature.^ There are two main 
changes. First, the degree of coUinearity needed between the quark and the photon in order 
for the pinching-pole approximation to be valid is modified. Since the thermal width is 
0{g'^T), one might expect that 6E ~ g'^T would be required for the poles to pinch. As 
6E ~ ^(pi/Pii^), smaller k means that p^ may be larger. If ~ g'^T, this would suggest 
that the poles still pinch for p^ ~ T^. Since we relied on p^^ <^ py^ at many points in our 
analysis, it appears that our treatment may break down when k < g^T. 

In fact this estimate is too pessimistic. As k becomes smaller, the cancellation we found 
which limits sensitivity to very-soft scattering events with |q| <C gT begins operating at 
larger values of |q|. In fact, there is an 0{q^/p]_) suppression for scatterings with <^ p^. 
The pinching pole requirement is that 6E be comparable to F, but with F given by the 
thermal width due to scatterings with momentum transfer larger than the characteristic 
|p_l| (not gT). Hence, the effective size of F is reduced to O {g'^ T^ / p'jj , up to logarithms. 
The typical value of p^^, for which F ~ 6E, is therefore p^ ~ g'^T'^^T /k. This means that 
our treatment only breaks down when the photon momentum reaches the scale g'^T, up to 
logarithms. This coincides with the inverse mean free path on which quarks undergo large 
angle scatterings. Extending our analysis below this scale would require major changes. 

One other change required for sufficiently small k is the correct inclusion of the modified 



•^^If the photon momentum is instead large compared to the temperature, fc 3> T, then the emission rate 
becomes exponentially small due to the statistical distribution functions, but all of our results remain valid. 
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photon dispersion relation which results from photon interactions with the plasma. Forward 
scattering off the plasma increases the photon frequency so that {k^—k) ~ e^T'^/k (or 
~ e^T^). This correction is non- negligible and must be included unless {k^—k) is much 
smaller than 5E for the typical range relevant to photon production. The analysis above 
shows that 5E ~ g'^y/kT, and hence corrections to the photon dispersion relation can be 
neglected only when e^T'^/k <^ g^y/kT, or k/T ^ Our results remain valid 

for photon momenta in the range eT <^ k< {e/gY/^T, provided one solves the integral 
equation ( |1.6|) with the correct photon dispersion correction included in 5E. If /c < eT, 
then plasma effects so strongly affect the propagation of the photon that it no longer travels 
at nearly the speed of light (and is more properly viewed as a plasmon — a charge density 
wave in the plasma). All of our analysis breaks down in this region. 
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•^^Note however that if we were considering gluon emission, then dispersion corrections due to the medium 
would always be important, even for k ^ T, a point which has been ignored in some previous treatments 
Q , although not in others Q . 
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